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Introduction 


It  is  now  well  established  that  self-assembled  (In,Ga)As  quantum  dots  grown  on  GaAs 
substrates  can  absorb  infrared  radiation  in  the  spectral  region  from  5.0  to  14  pm.  The 
absorption  is  due  to  inter-sublevel  transitions  in  the  conduction  band  of  the  (In,Ga)As 
structures.  When  used  in  the  active  regions  of  light  sources,  properly  designed  quantum 
dots  can  be  caused  to  emit  light  from  about  890  through  1350  nm.  The  access  to  these 
extended  spectral  regions  that  the  quantum  dots  provide,  especially  when  the 
technologically  advanced  and  robust  GaAs-based  materials  system  is  used,  has 
generated  a  lot  of  excitement  and  interest  in  the  area  of  quantum  dot  research.  A 
number  of  problems  still  remain  after  the  preliminary  demonstration  of  devices  based  on 
these  nanostructures. 

Statement  of  the  Problem 

The  potential  for  the  use  of  the  quantum  dots  in  mid-  infrared  detectors  and  near-infrared 
lasers  remains  to  be  realized.  Our  initial  results,  and  those  from  other  researchers  in  the 
field,  have  revealed  a  number  of  problems  associated  with  the  use  of  the  dots  in 
photodelectors  and  lasers.  The  key  problems  are  connected  with  (i)  uniformity  and 
control  of  the  dot  size,  (ii)  optimization  of  device  structures  and  the  spatial  location  of 
dopant  impurity  atoms,  and  (iii)  lack  of  in-plane  lateral  ordering  of  the  dots.  These 
problems  manifest  themselves  in  a  number  of  ways  in  real  devices.  In  lasers,  the  lack  of 
size  uniformity  leads  to  irreproducible  wavelengths  of  emission,  and  in  photodection,  it 
leads  to  broad  and  variable  absorption  lines.  The  lack  of  lateral  order  leads  to 
unpredictable  threshold  currents  in  lasers.  In  detectors,  dot  non-uniformity  and  improper 
doping  lead  to  detector  noise  and  hence  low  detectivity. 

Summary  of  Kev  Findings 

During  the  research  period  covered  in  this  report,  we  studied  the  structural 
characteristics  of  the  quantum  dots,  both  theoretically  and  experimentally.  The 
theoretical  analysis  was  undertaken  to  understand  the  important  physical  variables  (size, 
mole  fraction,  strain)  and  how  they  control  the  electronic  structure  and  hence  the 
emission  or  absorption  wavelengths  do  the  dots.  The  experimental  study  was  undertaken 
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to  validate  the  theory,  and  with  guidance  from  the  analytical  understanding,  find  out 
new  ways  to  better  engineer  and  build  the  structures  for  use  in  device  applications. 
Experimentally,  it  was  found  that  ordering  of  the  quantum  dot  nanostructures  could  be 
induced  in  the  vertical  direction  on  both  high  index  surfaces  as  was  as  on  the 
conventional  (001)  surface.  Lateral  ordering  was  much  more  difficult  on  either  substrate 
surface. 

The  first  key  result  is  from  our  theoretical  work  and  it  was  the  demonstration  that 
a  semi-analytical  method  could  be  used  to  provide  the  energy  structure  of  the  dots.  The 
method  we  developed  is  described  in  detail  in  Reference  [1].  A  copy  of  the  published 
paper  is  attached  to  this  report.  The  appeal  of  the  method  we  developed  is  that  solutions 
to  the  band  structure  of  the  dots  can  be  rapidly  obtained  without  recourse  to  the 
traditional  methods  that  require  intensive  computational  effort  [2,3],  including  in  some 
cases,  the  use  of  a  super  computer.  The  solutions  obtained  using  our  method  are  found 
to  be  in  compete  agreement  with  those  obtained  using  the  more  computationally  intensive 
methods.  Furthermore,  our  theoretical  results  agree  with  our  experimental 
measurements.  The  method  we  have  developed  therefore  has  become  a  useful  tool  in  our 
experimental  investigations  of  devices  fabricated  from  (In,Ga)As  quantum  dots.  We 
have  confirmed  both  experimentally  and  theoretically  that  most  dots  grown  in  our 
laboratory  are  pyramidal  in  shape,  but  with  “truncated”  tops. 

A  second  key  result  was  from  our  experimental  work.  We  demonstrate  that  by 
embedding  quantum  dots  in  symmetrically  graded  quantum  wells,  the  near-infrared 
emission  of  the  dots  can  be  enhanced  by  as  much  as  an  order  of  magnitude.  We  further 
showed  that  the  emission  wavelengths  of  such  dots  can  be  extended  to  the  1300  nm 
range,  which  is  a  of  technological  importance.  Silica-based  fibers  for  communication 
exhibit  an  absorption  and  a  dispersion  minimum  at  1300  nm.  The  results  have  been 
reported  n  the  literature  [4]  and  a  copy  of  the  published  paper  is  attached  to  this  report. 

A  third  key  finding  from  our  work  is  that  when  the  (In,Ga)As  quantum  dots  are 
used  to  fabricated  infrared  inter-sublevel  detectors,  the  noise  in  these  devices  is  primarily 
due  to  generation  -recombination  noise  between  78  K  -  100  K.  This  can  be  reduced  by 
improving  the  quality  of  the  structures.  In  particular,  one  can  obtain  results  by 
improving  the  electron  mobility  and  lowering  the  defects  in  die  device  structures. 
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The  fourth  finding  is  that  the  quantum  size  variation  leads  to  broad  and  random 
absorption  lines.  These  results  have  been  reported  in  the  technical  literature  [5]  and  a 
copy  of  the  published  paper  is  attached  to  this  report. 
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Abstract 

This  paper  presents  a  simple  analytical  method  for  calculating  the  strain  distribution  in  and  around  self-assembled 
(In,Ga)As/GaAs  quantum-dot  nanostructures.  The  dots  are  assumed  to  be  buried  in  an  infinite  medium  so  that  the  effects  of 
tree  surfaces  can  be  neglected.  This  assumption  is  based  on  the  relative  size  of  the  dot,  compared  to  that  of  the  overlayer. 
The  model — based  on  classical  continuum  elasticity — is  capable  of  handling  dots  of  arbitrary  shapes;  here,  however,  only 
dots  with  pyramidal  and  truncated-pyramidal  shapes  are  considered.  The  approximate  shape  of  the  dots  is  extracted  from 
high-resolution  transmission  electron  microscope  observations.  The  electronic  energy  levels  in  the  dots  are  calculated  by 
solving  the  three-dimensional  effective  mass  Schrodinger  equation.  The  carrier  confinement  potential  in  this  equation  is 
modified  by  the  strain  distribution.  Because  the  dots  are  in  a  strong  confinement  regime,  the  effects  of  Coulomb  interactions 
are  neglected.  The  calculated  confined  eigen-energies  agree  with  our  experimental  photoluminescence  data.  The  calculations 
also  support  previous  results  reported  by  others.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 

PACS:  83.85. St;  73.21. La;  68.66.Hb;  68.37.Lp;  73.43.Cd 

Keywords:  Quantum  dots;  Strain;  Eight-band  model;  Electronic  spectra;  Nanostructures 


1.  Introduction 

Progress  in  nanotechnology  has  resulted  in  the 
synthesis  of  nanometer-scale  self-assembled  quantum 
dots  for  potential  applications  in  electronics  and  op¬ 
toelectronics.  Several  devices  that  include  quantum 
dots  in  their  active  regions  have  been  proposed  in 
the  past  few  years;  some  of  them,  such  as  lasers, 
infrared  detectors,  and  memory  elements  have  been 
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demonstrated.  Since  millions  of  dots  are  generally 
required  to  yield  macroscopic  effects,  a  narrow-size 
distribution  is  important  for  optimum  device  perfor¬ 
mance.  The  synthesis  of  uniform-sized  dots  is  still 
a  research  problem  and  much  remains  to  be  learned 
in  understanding  the  growth  processes  that  lead  to 
pristine  and  uniform-sized  quantum  dots.  The  most 
studied  quantum-dot  system  today  is  based  on  the 
spontaneous  formation  of  (In,Ga)As  islands  on  a 
GaAs  substrate  during  heteroepitaxy.  This  process  is 
a  consequence  of  the  strain  which  results  from  the 
in-plane  lattice-mismatch  between  the  (In,Ga)As  and 
the  GaAs.  The  specific  shape  of  the  dots  is  usually 


1386-9477/02/$ -see  front  matter  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
PII:  SI  386-9477(02)00459-9 


132 


V-G.  Stoleru  et  al  / Physica  E  15  (2002)  131-152 


a  function  of  the  growth  parameters  and  the  general 
environment  of  the  dots.  Some  researchers  have  ob¬ 
served  dots  with  hemispherical  shapes;  others  have 
observed  dot-shapes  that  are  multifaceted  domes; 
and  yet  several  other  groups  have  reported  dots  with 
pyramidal  shapes  [1].  It  is  generally  difficult  to  de¬ 
termine  the  shape  of  the  dots  in  situ  during  growth; 
the  lateral  extent  and  height  of  the  dots,  however, 
can  often  be  determined  approximately  for  surface 
quantum  dots  from  atomic  force  microscope  stud¬ 
ies.  Better  estimates  of  buried  dot-sizes  can  usually 
be  obtained  from  transmission  electron  microscope 
studies;  these  studies  can  generally  provide  valuable 
information  about  the  shape  of  the  dots  [2].  Central  to 
understanding  many  properties  of  the  dot  is  a  need  to 
determine  the  elastic  strain  distribution  in  and  around 
the  dots.  The  strain  distribution  profoundly  affects 
the  electronic  structure,  and  hence  the  optical  proper¬ 
ties  of  the  dots  by  modifying  the  energies  and  wave 
functions  for  confined  carriers.  The  principal  aim  of 
this  paper  is  to  quantitatively  determine  the  stress 
fields  and  hence  the  strain  distribution  in  and  around 
a  dot.  This  information  is  then  used  to  calculate  the 
electronic  structure  of  the  dots. 

Strain  drives  the  Stranski-Krastanow  growth  pro¬ 
cess  that  leads  to  the  formation  of  dots;  residual  strain, 
in  general,  can  have  interesting  effects  on  the  elec¬ 
tronic  and  optical  properties  of  semiconductor  struc¬ 
tures.  The  effects  of  strain  on  semiconductors  have 
therefore  motivated  a  significant  interest  in  calculating 
its  magnitudes  in  a  broad  range  of  devices  [3].  The  hy¬ 
drostatic  (Sh)  component  of  strain,  for  example,  usu¬ 
ally  shifts  the  conduction  and  valence  band-edges  of 
semiconductors;  biaxial  (eb)  strain,  on  the  other  hand, 
modifies  the  valence  bands  by  splitting  the  degener¬ 
acy  of  the  light-  and  heavy-hole  bands.  These  effects 
have  a  profound  impact  on  the  electronic  and  optical 
properties  of  the  structures  out  of  which  devices  might 
be  fabricated. 

The  two  general  classes  of  techniques  used  in 
the  calculation  of  strain  include  the  finite  difference 
[4,5]  and  atomistic  calculation  methods  [6].  Both  of 
these  methods  require  considerable  computational 
resources.  In  the  atomistic  approach,  there  is  an  im¬ 
plicit  reliance  on  the  validity  of  the  valence  force 
field  model  with  a  Keating  interatomic  potential. 
This  method,  however,  has  recently  been  shown  to 
be  questionable  for  small  islands  whose  base  dimen¬ 


sions  are  smaller  than  lOnm  [7];  for  these  islands, 
the  Stillinger-Weber  potential  has  been  suggested  as 
an  alternative  potential.  The  finite  element  and  finite 
difference  methods  have  been  used  by  several  groups 
[4,8-11].  One  advantage  of  the  finite  element  and 
finite  difference  methods  is  that  they  can  be  used  to 
treat  problems  with  complex  geometries.  This  advan¬ 
tage,  however,  is  offset  by  the  demand  on  computa¬ 
tional  resources.  Another  variant  of  these  methods  is 
the  boundary  element  approach;  this,  however,  can  be 
mathematically  complex  [12].  A  simple  and  elegant 
method  for  calculating  strain  fields  around  a  single, 
isotropic,  cubic  dot  has  been  presented  by  Downes 
et  al.  [13].  This  method  is  based  on  a  simplification  of 
Eshelby’s  classic  inclusion  theory  [14].  The  method 
first  identifies  a  set  of  vectors  such  that  the  divergence 
of  each  gives  the  Green’s  function  for  the  stress  com¬ 
ponents  ay.  By  invoking  Gauss’  theorem,  the  stress 
field  is  determined  by  performing  an  integral  over  the 
surface  of  the  dot. 

In  this  paper,  we  extend  the  method  to  more  com¬ 
plex  and  practical  geometries.  We  particularly  focus 
on  regular  pyramids  and  those  with  truncated  tops. 
In  the  model,  we  make  the  basic  assumption  that  the 
elastic  properties  of  the  materials  are  isotropic  [15]. 
This  assumption  does  not  seriously  affect  any  qualita¬ 
tive  conclusions  drawn  from  the  calculations.  In  any 
event,  it  has  been  shown  in  a  recent  paper  [16]  that  as 
long  as  the  symmetry  of  the  shape  of  the  structure  is 
less  than  or  equal  to  the  cubic  symmetry  of  the  crys¬ 
tal,  both  anisotropic  and  isotropic  models  give  similar 
results.  We  caution  however  that  even  though  this  is 
true  for  quantum-dot  nanostructures,  it  is  not  so  for 
quantum- well  heterostructures. 

Our  calculations  here  are  for  the  (In,Ga)As/GaAs 
materials  system.  The  model,  however,  is  general 
enough  that  it  can  be  used  for  similar  other  materi¬ 
als  systems.  For  the  InAs/GaAs  system,  the  initial 
strain  or  lattice-mismatch  is  defined  as  £o  =  (^GaAs  — 
0inAs)MnAs  =  —0.067.  This  is  the  fractional  change  of 
the  difference  between  the  in-plane  lattice  constant  of 
the  GaAs  substrate  and  the  InAs  epilayer  with  respect 
to  the  InAs  lattice  constant;  when  Ga  is  included  in 
the  composition  of  the  epilayer,  the  mismatch  (strain) 
becomes  £0  ”  O^GaAs  -  ^inpai_^As)/ain;(Gai_xAs. 
The  Poisson  ratios  for  the  binary  materials  are  taken 
to  be  0.316  for  the  GaAs  and  0.354  for  the  InAs.  As 
is  usual  in  continuum  elasticity,  we  assume  that  the 
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materials  are  continuous,  linear,  isotropic,  and  obey 
Hooke’s  law.  The  use  of  continuum  elasticity  meth¬ 
ods  has  been  experimentally  verified  to  be  valid  for 
layers  as  thin  as  three  atomic  monolayers  [17];  fur¬ 
thermore,  it  is  estimated  theoretically  that  the  method 
remains  valid  for  layers  whose  thicknesses  are  on 
the  order  of  five  atoms  [18].  We  want  to  point  out 
that  the  sensitivity  of  some  physical  properties  to 
strain  could  make  anisotropic  effects  important;  the 
isotropic  approximation  should  therefore  be  treated 
with  caution,  particularly  for  layers  oriented  in  certain 
crystallographic  directions.  For  most  cases,  however, 
the  anisotropy  only  modifies  the  strain  distributions 
slightly  [16,19]. 

Classical  elasticity  and  atomistic  methods  have 
been  shown  to  give  similar  results  for  small  strains 
[20].  Differences  become  apparent  for  strains  larger 
than  5%,  particularly  in  the  case  of  semiconduc¬ 
tor  quantum  dots.  Even  though  atomistic  elasticity 
is  expected  to  continue  to  give  reliable  results  for 
large  strains,  the  number  of  atoms  involved  in  the 
computations  makes  the  method  unwieldy  for  rapid, 
order-of-magnitude  calculations.  The  advantage  of 
the  approach  based  on  Eshelby’s  inclusion  theory  is 
that  it  is  simple  and  allows  one  to  find  approximate 
analytical  expressions  for  the  strain  tensor  compo¬ 
nents.  As  we  will  show  later,  for  quantum  dots  with 
lateral  dimensions  in  the  range  of  10-20  nm,  the 
strain  profiles  calculated  within  the  modified  Eshelby 
framework  are  in  good  agreement  with  those  obtained 
by  the  methods  of  atomistic  elasticity. 

This  paper  is  organized  as  follows:  in  Section  2 
we  give  a  detailed  description  of  the  basic  model  for 
calculating  the  strain  tensor  for  pyramidal  and  trun¬ 
cated  pyramidal  quantum  dots.  The  strain  fields  are 
presented  and  discussed  in  Section  3.  The  method  dis¬ 
cussed  here  is  a  suitable  precursor  for  the  calculation 
of  electronic  structure  based  on  the  envelope  function 
method  within  the  plane-wave  expansion  technique. 
One  can  determine,  for  example,  the  strain-dependent 
matrix  element  linking  any  pair  of  plane  waves.  We 
present  and  discuss  the  strain-dependent  confining  po¬ 
tentials  for  electrons  and  holes,  as  well  as  the  elec¬ 
tronic  energy  levels  in  Section  4.  A  summary  of  our 
results  is  given  in  Section  5. 

In  our  calculations,  we  consider  a  single  isolated  dot 
to  facilitate  comparison  with  published  results.  The  dot 
is  assumed  to  be  buried  deep  within  a  matrix  material. 


Figs.  1(a)  and  (b)  show  the  schematic  cross-sections 
of  the  pyramidal  and  truncated  pyramidal  quantum-dot 
structures  under  consideration.  They  are  assumed  to 
be  on  top  of  a  semi-infinite  (001)  GaAs  substrate  on 
which  is  grown  a  thin  InAs  wetting  layer  first. 


2.  Model  description 


Most  theoretical  calculations  of  the  properties  of 
InAs  dots  assume  a  square-based  pyramidal  shape 
[4,6,21].  Here,  for  purposes  of  comparison  with  previ¬ 
ously  published  work,  we  also  consider  this  geometry. 
However,  we  extend  our  calculations  to  include  the 
truncated  pyramidal  case.  When  a  capping  layer — for 
example  GaAs — is  grown  on  top  of  a  layer  of  dots, 
the  morphology  of  the  overlayer  is  affected  by  the  in¬ 
teraction  of  the  inhomogeneous  strain  around  the  dots 
and  in  the  wetting  layer.  This  affects  the  apex  of  the 
dots.  There  is  then  a  thermodynamically  favored  ten¬ 
dency  for  the  adatoms  to  migrate  to  the  side  of  the 
dots  [22],  resulting  in  a  reduction  of  the  dot  height  and 
formation  of  a  flat  (0  0  1)  top  surface. 

Following  Downes  et  al.  [13],  the  Lame  potential  u 
during  relaxation  can  be  described  by  a  scalar  potential 

(1) 

where  G  is  the  shear  modulus,  defined  as  2G~E/(\  + 
v)  in  Ref.  [23].  The  displacement  potential  obeys  the 
Poisson  equation 

V2$(r)  =  So(r)2G.  (2) 


From  Eq.  (2),  a  solution  can  be  reached  by  integration, 
using  Green’s  function,  thus 


(1+v) 

(1  -  V) 


gpOo) 
47t|r  -  ro| 


d  3r0. 


(3) 


In  the  equation  above,  the  function  l/|r  —  ro|  can  be 
written  as 


1^  (r-r q) 

2  |  r-  r0 1  ‘ 

The  point  /*o  is  within  the  volume  of  the  dot.  The  vol¬ 
ume  integral  in  Eq.  (3)  can  be  converted  to  a  surface 
integral  by  applying  the  divergence  theorem  to  it,  with 
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B=2a 

(c) 

Fig.  1.  Schematic  models  for  the  pyramidal  (a)  and  truncated  pyramidal  (b)  quantum-dot  structures;  2?  =  base  width,  H  —  pyramid  height, 
/  =  truncation  factor,  h  —  truncated  pyramid  height,  d  —  thickness  of  wetting  layer;  (c)  Three-dimensional  model  of  a  quantum  dot  in  the 
Cartesian  coordinate  system. 


the  point  ro  now  being  on  the  surface  of  the  dot,  as 
shown  in  Fig.  1(c).  Thus 


—  <p(r)  =  —  OJlll 

2 G  K  ’  in  (1  -  v) 


\r  -  ro\ 


1  (1  +  v) 

8rc  (1  -  v) 

* II 


(4) 


The  initial  lattice-mismatch,  eo>  can  be  considered 
constant  within  the  volume  of  the  dot.  With  the  ap¬ 
propriate  substitutions,  Eq.  (1)  becomes 


u=j -V0=±(i±i>  if 

2 G  4tt  (1  -v)JJs 


\r-r0 


dS(r<)).  (5) 


The  stress  components  can  therefore  be  written  as 
£qE 


Vij(r)  ■■ 


87t(1  —  v) 

[  [  (i-  io)i  +  (j  ~  joV  4  r/„  s 

Jl — F7^!3 — <W('*) 


+  <5; 


(6) 


where  i  and  j  are  unit  vectors  in  the  zth  and  yth  di¬ 
rections,  respectively,  for  i  =  x,  y,  z,  j  =  x,y,z  and 
r2=x 2+y2jt-z2.  The  parameter  so  is  the  isotropic  misfit 
strain,  E  is  the  Young’s  modulus,  v  is  Poisson’s  ratio, 
and  Sjj  is  the  Kronecker  delta  function.  The  last  part 
of  Eq.  (6)  comes  from  evaluating  the  limit  of  the  sur¬ 
face  integral,  as  a  field  point  r  approaches  the  bound¬ 
ary  point  ro  on  the  surface  of  the  dot.  The  misfit  strain 
is  taken  as  negative  for  a  material  under  compression. 
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The  volume  of  a  square-based,  truncated  pyramid  is 
defined  by 


-B 


2 H 
(H  —  z) 


2 H 
XH-z) 


<  y  <5- - -, 

2 H  7  2H 

0  <  z  <  Ht, 


(7) 


and 


a'zz  =  4n  —  4 


—tan 


-l 


z  —  a 
\z  -  a\ 


+  tan 


\/2a2  +  z2/ J 


+ 


In 


VTflT+z2  —  a 


\/2a2  4-  z2  +  a 


-In 


\/2fl2  +  z2  +  a 


V2a2  -I-  : 


where  H  is  the  height  of  the  pyramid  in  the  absence  of 
truncation,  B  is  the  base  of  the  pyramid,  and  0  t  <  1 , 
where  t  represents  the  truncation  factor.  The  z-axis 
is  the  [001]  growth  direction,  and  the  origin  of  the 
coordinates  is  at  the  center  of  the  square  base  of  the 
pyramid  (the  z  =  0  plane). 

After  converting  to  Cartesian  coordinates,  the  in¬ 
tegrations  indicated  in  Eq.  (6)  can  be  easily  carried 
out.  The  in-plane  stress  components  have  cumbersome 
mathematical  expressions;  however,  for  some  parti¬ 
cular  directions,  they  can  be  simplified.  For  example, 
for  a  square-based  pyramidal  dot,  with  a  contrast  ratio 
(defined  as  the  ratio  of  the  height  to  the  base)  of  1:2, 
where  B  —  2a  and  H  =  a,  one  obtains  the  following 
expressions  as  functions  of  z: 


Note  that  a,y  =Eso/(4n(l  — v))cjA,  where  So,  v,  and  is 
are  as  defined  before.  The  more  general  analytical  ex¬ 
pressions  for  calculating  stress  distributions  in  struc¬ 
tures  with  arbitrary  degrees  of  truncation  are  given  in 
the  Appendix  A.  The  expressions  for  the  strain  com¬ 
ponents  follow  by  substitution  of  the  stress  compo¬ 
nents  into  Hooke’s  law.  The  generalized  stress-strain 
relations  are: 

Vkl  =  CklmnZmn,  (10) 


^XX  G  yy 


=  2 
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x/2  a2  +  z2  —  a 
(z  -  a) 


In 


v/3 


+  \z-a\ 


/— 2  a  —  z 

ini— - + 


\f 2  a2  +  z2^ 


(8) 


and 

£kl  =  SklmnGrr, 


(11) 


where  the  Q/„;„  are  the  compliance  and  the  Sum„  the 
stiffness  coefficients.  For  materials  with  cubic  symme¬ 
try,  only  three  of  the  81  components  are  independent. 
So  Ci i  =  Ci in  =  C2222  —  C3333,  C\2  -  C1122  =  C2233, 
etc.,  and  C44  =  C1212  =  C2323,  etc.;  similar  relations 
hold  for  Sn,  1S12  and  .S'44.  The  stress-strain  relations 
can  then  be  written  in  terms  of  Young’s  modulus  E 
and  Poisson  ratio  v.  Since 


E  = 


1  (C„-Ci2XCn+2Ci2) 


Sn 


and 

V  =  -ES\2  = 


C\ 1  +  C12 

C12 


Cn+Cn' 


then 


■  %  + 


vE 


1+v  J  (1  +  v)(l  —  2v) 


&ij&kk: 


(12) 


(13) 


(14) 
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and 

fiy  =  ^[(l+v)ffy-  SijVffu].  (15) 

In  this  paper  we  define  the  hydrostatic  and  the  biaxial 
strains  as 


£h  —  &XX  +  £yy  4"  £zz> 

(16) 

and 

£b  =  &zz  2  (f**  ® yy )> 

(17) 

respectively.  These  relations  are  introduced  here  be¬ 
cause  they  represent  important  quantities  used  in  the 
analysis  of  the  electronic  energy  levels  in  a  later 
section. 


3.  Strain  distribution 

We  have  calculated  the  strain  distributions  for  a 
number  of  structures  whose  dimensions  are  given  in 
Table  1.  Calculations  for  pyramidal  QDs,  labeled  PQD 
in  Table  1,  as  well  as  truncated  pyramidal  dots,  labeled 
TPQD  in  Table  1,  have  been  carried  out.  We  consider 
first  structure  PQD3  in  Table  1;  this  type  of  structure 
has  been  extensively  studied  and  reported  on  in  the 
literature.  Fig.  2  shows  the  strain  tensor  components 
Sxx  and  £zz  for  structure  PQD3,  plotted  as  functions 
of  position  along  the  z-axis.  The  shear  strain  compo¬ 
nents,  exy ,  Syzi  and  £*z  are  negligible  in  the  dot  and 
barrier  materials;  they  could,  however,  be  apprecia¬ 
ble  at  the  interfaces  [4].  By  symmetry,  the  component 
Bxx  =  eyy.  In  a  thin  substrate  region  below  the  dot, 


the  GaAs  lattice  experiences  a  tensile  (positive  £**) 
strain  in  the  x-y  plane  and  a  compressive  (negative 
£zz)  strain  in  the  z  direction.  In  this  case,  the  dot  is 
forcing  the  substrate  lattice  constant  to  be  that  of  InAs 
(tfinAs  —  6.05  A).  In  the  base  region  of  the  dot,  the  situ¬ 
ation  is  reversed.  Here,  ezz  is  positive  and  £**  negative 
because  the  substrate  now  attempts  to  force  the  dot 
lattice  constant  to  be  that  of  GaAs  (tfQaAs  =  5.65  A). 
With  increasing  height  within  the  dot,  szz  changes  its 
sign,  becoming  negative  at  the  top  of  the  pyramid.  At 
the  top  of  the  pyramid,  the  dominant  forces  acting  on 
the  dot  originate  from  the  GaAs  matrix  at  the  sides, 
causing  a  compressive  strain  (negative  ezz)  along  the 
z  direction  and  a  tensile  strain  (positive  in  the 
x-y  plane.  Fig.  3  is  a  plot  of  the  e**  and  szz  components 
along  the  x-axis  for  structure  PQD3.  Within  the  dot, 
both  £xx  and  are  negative,  implying  the  existence 
of  a  region  of  hydrostatic  compression.  The  £**,  eyy , 
and  £zz  components  of  the  strain  tensor  are  plotted  in 
Figs.  4(a),  (b),  and  (c)  in  the  x-z  plane. 

The  hydrostatic  and  biaxial  components  of  the  strain 
for  structure  PQD3  are  plotted  as  functions  of  posi¬ 
tion  along  the  z-axis  in  Fig.  5.  The  hydrostatic  strain 
is  compressive  within  the  dot  and  in  the  surrounding 
barrier  material.  The  interior  of  the  quantum  dot  ex¬ 
hibits  a  nearly  homogeneous  hydrostatic  strain,  while 
the  barrier  experiences  a  small  hydrostatic  strain.  This 
is  because  GaAs  is  stiffer  than  InAs.  The  biaxial  strain 
tends  to  be  negative  in  the  barrier  and  positive  in  the 
dot,  and  it  is  zero  near  the  center  of  the  dot.  In  this 
region,  the  strain  is  entirely  hydrostatic  in  character. 
Fig.  5  suggests  that  a  significant  transfer  of  biaxial 
strain  to  the  barrier  takes  place. 


Table  1 

Dimensions  of  single  (In,Ga)As/GaAs  quantum-dot  structures  used  for  computations 


Structure 

Base  width 

B  (nm) 

Pyramid  height 

H  (nm) 

Truncation 
factor  t 

Truncated  pyramid 
height  h  (nm) 

PQD1 

12 

3 

PQD2 

12 

4 

PQD3 

12 

6 

PQD4 

16 

6 

PQD5 

20 

6 

TPQD1 

13.8 

4.6 

0.25 

1.15 

TPQD2 

13.8 

4.6 

0.50 

2.30 

TPQD3 

13.8 

4.6 

0.75 

3.45 

TPQD4 

19.7 

4.7 

0.64 

3.00 
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Fig.  2.  Strain  tensor  components  txx  and  ezz  for  structure  PQD3 
plotted  along  the  z-axis. 
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Fig.  3.  Strain  tensor  components  Exx  and  ezz  for  structure  PQD3 
plotted  along  the  jc-axis. 


We  have  studied  strain  relaxation  in  the  pyramidal 
InAs/GaAs  QDs  for  the  different  structures  labeled 
PQD(l-5)  in  Table  1 .  We  compare  the  strain  distribu¬ 
tions  in  the  x-z  plane  for  dots  of  different  base  widths 
(B\  =  12  nm  for  PQD3,  2?2  =  16  nm  for  PQD4,  and 
£3 =20  nm  for  PQD5)  but  the  same  height  (H= 6  nm) 
in  Fig.  6.  We  note  that  for  dots  with  small  base  widths, 
the  strain  components  change  sign  from  the  base  re¬ 
gion  of  the  dot  to  the  apex  more  rapidly  than  for 
dots  with  large  base  widths.  For  example,  in  structure 
PQD5,  the  strain  components  £**  and  eyy  are  negative 
inside  the  dot,  while  in  structure  PQD3,  they  change 
sign  from  negative  at  the  base  of  the  dot  to  positive  at 
the  apex  of  the  dot.  The  strain  component  ezz  changes 
rapidly  from  positive  to  negative  values  along  the 


Fig.  4.  Strain  tensor  components  for  structure  PQD3  i 
plane:  (a)  Exx ,  (b)  £yy,  (c)  ezz  (strain  expressed  in  %). 
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Fig.  5.  The  hydrostatic  £h,  and  biaxial  components  of  the  strain 
for  structure  PQD3  plotted  along  the  z-axis. 


z-axis  in  structure  PQD3,  but  is  mostly  tensile  inside 
structure  PQD5.  The  electronic  structure  of  the  dots 
is  also  expected  to  have  such  a  strong  dependence  on 
the  base  lengths,  as  we  will  discuss  in  Section  4. 

Fig.  7  illustrates  the  strain  relaxation  pattern  for 
InAs  dots  of  the  same  base  width  (B  =  12  nm)  but 
different  heights  (H\  =  3  nm  for  PQD1,  i72  =  4  nm 
for  PQD2,  and  H$  =  6  nm  for  PQD3).  This  compari¬ 
son  suggests  that  the  strain-modified  electronic  energy 
structure  is  likely  not  to  be  as  sensitive  to  variations 
in  the  height  of  the  dot  as  it  is  to  changes  in  the  base 
width  of  the  dot.  Extraordinarily  tall  dots  have  large 
hydrostatic  strains  near  their  centers.  This  increases 
the  band  gap  inside  the  islands,  partially  compensat¬ 
ing  for  the  change  in  confinement  potential. 

Calculations  indicate  that  when  InAs  self-assembled 
pyramidal  dots  are  buried  under  a  GaAs  overlayer, 
significant  tensile  stress  is  induced  at  the  top  of  the 
dot;  the  stress  at  the  edges  of  the  dot,  on  the  other 
hand,  is  weak  and  compressive.  This  is  in  contrast 
to  surface  dots,  where  stress  relaxation  occurs  at  the 
island  tops  but  is  concentrated  at  the  edges  [10].  It 
should  be  noted  that  the  stress  component  (at  the 
x  =  0  plane)  is  compressive  (negative)  in  embed¬ 
ded  dots  whereas  it  is  tensile  (positive)  in  uncovered 
ones.  The  strain  distribution  maps  calculated  here  for 
pyramidal  InAs  dots  are  in  reasonably  good  agree¬ 
ment  with  published  results  obtained  by  calculations 
based  on  atomistic  methods  [4,24].  The  agreement 
is  surprisingly  good  considering  the  computational 
simplifications  introduced  in  the  present  work. 

We  have  extended  the  analytical  procedure  to  calcu¬ 
lating  strain  distributions  of  truncated  InAs  pyramidal 


dots  such  as  that  shown  in  the  transmission  electron 
microscope  (TEM)  image  of  Fig.  8.  The  structure  in 
this  figure  was  grown  by  molecular  beam  epitaxy  on 
a  (0  0 1 )  GaAs  substrate.  From  the  substrate  up,  the 
structure  consists  nominally  of  a  250-nm-GaAs  buffer 
layer;  this  is  followed  by  five  periods  of  InAs/GaAs 
quantum  dots.  The  entire  structure  is  caped  with  27  nm 
of  GaAs.  The  dot  density,  as  determined  by  atomic 
force  microscope,  was  about  4  x  1010  cm~2.  The  shape 
of  the  dots  was  determined  from  cross-sectional  TEM 
studies.  The  TEM  image  clearly  shows  that  the  dots 
are  coherent,  with  no  observable  dislocations.  The  dots 
could  be  discerned  from  their  darker  appearance  due 
to  the  presence  of  Indium.  The  buried  dots  appear  to 
be  pyramids  with  truncated  tops;  their  base  and  height 
were  found  to  be  B  =  13.8  nm,  and  h  =  3.4  nm,  re¬ 
spectively.  This  corresponds  to  a  truncation  factor  of 
t  =  0.75  (structure  TPQD3  in  Table  1).  This  observa¬ 
tion  is  consistent  with  previous  reports  [25].  We  have 
calculated  the  strain  tensor  components,  as  well  as  the 
hydrostatic  and  the  biaxial  strain  profiles  for  this  struc¬ 
ture.  The  results  are  shown  in  solid  lines  in  Fig.  9.  We 
mention  that  the  analytical  results  obtained  for  struc¬ 
ture  TPQD3  are  in  very  good  agreement  with  those 
obtained  for  a  similar  structure,  using  finite  element 
method,  in  Ref.  [26]. 

For  pyramidal  dots  of  arbitrary  truncation,  general 
expressions  for  the  stress  tensor  components  can  be 
derived.  These  expressions  are  given  in  the  Appendix 
A  for  a  uniformly  lattice-mismatched  InAs/GaAs  dot. 
Calculations  based  on  these  expressions  for  truncated 
InAs/GaAs  dots  have  been  performed  for  structures 
TPQD(l-3)  in  Table  1.  The  specific  truncation  fac¬ 
tors  used  are  t  =  0.25,  0.50,  and  0.75,  respectively. 
For  t  =  1  we  obtain  the  whole  pyramidal  geometry, 
while  for  t  =  0,  we  merely  reproduce  results  for  the 
two-dimensional  wetting  layer  alone. 

There  are  several  features  that  are  common  to  all 
strain  distributions  for  the  structures  studied.  First,  the 
magnitudes  of  the  strain  components  are  largest  at 
the  dot/matrix  interface,  particularly  at  the  vertices. 
The  strain  is  rapidly  attenuated  within  the  dot  and 
in  the  matrix  material  immediately  below  the  square 
base  of  the  dot.  It  remains  of  similar  form  for  each 
truncation  since  the  base  of  the  dot  is  unchanged.  The 
hydrostatic  strain  is  zero  in  the  matrix  material  and 
proportional  to  So  in  the  quantum  dot,  consistent  with 
standard  inclusion  theory.  The  conduction  band  shift 


z[nm]  z[nm]  z[nm] 


V.-G.  Stoleru  et  al.  /  Physica  E  15  (2002)  131-152 


Fig.  6.  Comparison  of  strain  components  for  structures  PQD(3-5),  characterized  by  height  H  =  6  nm,  but  different  base  widths:  B\  —  12  nm 
(PQD3),  B2  =  16  nm  (PQD4),  and  B3  —  20  nm  (PQD5):  (a)  exc,  (b)  eyy ,  (c)  Ezz  (strain  expressed  in  %). 
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Fig.  6.  ( Continued. ) 


depends  on  the  hydrostatic  strain,  so  the  conduction 
band  shift  is  relatively  constant  within  a  dot  of  constant 
composition.  Figs.  9(a)  and  (b)  illustrate  the  sxx  and 
ezz  strain  tensor  components  plotted  along  the  z-axis 
for  structures  TPQD(1~3)  of  Table  1. 

Our  results,  and  those  of  others,  indicate  that  the 
magnitudes  of  the  strain  components  depend  on  the 
geometries  of  the  dots.  This  is  as  it  should  be  for 
both  qualitative  and  quantitative  reasons.  Consider, 
for  example,  that  the  biaxial  strain  is  very  sensitive  to 
truncation  (see  Fig.  9(c));  at  the  apex  of  the  pyramid, 
the  biaxial  strain  is  negative,  but  becomes  increasingly 
positive  as  the  degree  of  truncation  increases. 

For  structure  TPQD4  we  consider  the  case  where 
the  dot  material  is  a  ternary  compound;  in  particular, 
we  consider  the  case  where  the  Indium  composition 
is  23%,  as  determined  from  X-ray  diffraction  experi¬ 
ments.  And  as  before,  the  shape  and  size  of  the  dots 
are  determined  from  cross-sectional  high-resolution 
TEM  observations.  We  have  performed  calcula¬ 
tions  to  determine  the  strain  tensor  components  for 
Ino.23Gao.77As/GaAs  dots.  The  strain  distributions  are 
depicted  in  Fig.  10  for  structure  TPQD4.  The  strain 
tensor  components  obtained  here  are  used  as  input  for 
the  electronic  band  structure  calculations  discussed  in 
the  next  section. 


4.  Electronic  structure  and  optical  transitions 

The  band  structure  of  semiconductors  is  generally 
altered  by  the  presence  of  strain,  which  changes  the 
lattice  constant  and  reduces  the  symmetry  of  the  crys¬ 
tal.  Strain  modifies  energy  gaps  and  removes  degener¬ 
acy.  Here,  we  take  account  of  the  strain  distributions 
discussed  in  the  previous  sections  in  calculating  the 
electronic  structure  of  dots  in  the  envelope  function 
approximation  using  an  eight-band  strain-dependent 
k  •  p  Hamiltonian  [27].  The  eight-band  k  •  p  method 
represents  an  extension  of  the  Luttinger-Kohn  formal¬ 
ism,  which  describes  coupling  among  the  light-hole, 
heavy-hole  and  split-off  valence  bands  to  second  order 
in  k,  but  is  modified  to  include  the  linear  coupling  be¬ 
tween  the  conduction  and  valence  band  states.  This  is 
necessary  in  order  to  correctly  model  conduction  band 
non-parabolicity.  A  product  of  strain  components 
and  a  deformation  potential  describes  the  general 
form  of  the  strain-induced  modification  of  the  band 
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Fig.  7.  Comparison  of  strain  components  for  structures  PQD(l-3),  characterized  by  base  width  £=12  nm,  but  different  heights:  H\  =3  nm 
(PQD1),  H2  —  4  nm  (PQD2),  and  Hi -6  nm  (PQD3):  (a)  e**,  (b)  eyyi  (c)  ezz  (strain  expressed  in  %). 
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Fig.  7.  ( Continued. ) 


Fig.  8.  Cross-sectional  high-resolution  transmission  electron  mi¬ 
croscope  image  of  InAs  quantum  dots  in  a  GaAs  matrix. 


structure.  Typical  values  of  deformation  potentials 
range  from  1  to  10  eV  [22].  These  determine  the 
amount  of  band  structure  modification  and  enter  the 
calculations  as  material  parameters  given  in  Table  2. 
In  our  calculations,  we  neglect  the  lack  of  inversion 
symmetry  in  the  zincblende  structure  and  conse¬ 
quently  ignore  the  small  coupling  of  the  conduction 
band  to  shear  deformations.  In  the  strain-dependent 
Hamiltonian,  the  non-diagonal  terms  containing  shear 
components  of  the  strain  are  non-zero  in  our  geome¬ 
tries,  except  far  away  from  the  dot  axis.  However, 
since  the  diagonal  shear  strain  term  of  the  Hamiltonian 
decouples  the  heavy-  and  light-hole  bands  by  approx¬ 
imately  100  meV,  the  correction  of  the  non-diagonal 
terms  in  the  Hamiltonian  is  small  near  the  band-edge. 
The  resultant  Hamiltonian  is  a  matrix  that  is  diag¬ 
onalized  using  the  Lanczos  algorithm  developed  in 
Refs.  [28,29], 

Some  insight  can  be  gleaned  by  examining  the 
strain-induced  modification  to  the  band  structure 
corresponding  to  structure  PQD3  in  Table  1.  In  the 
absence  of  strain,  the  confining  potential  for  an  elec¬ 
tron  (hole)  is  a  square  well  formed  by  the  difference 
in  the  absolute  energy  of  the  conduction  (valence) 
band-edges  in  InAs  and  GaAs  [30].  In  the  conduc¬ 
tion  band,  the  depth  of  the  confining  potential  in  this 
case  is  about  840  meV.  For  holes,  the  potential  well 
is  about  263  meV  deep.  However,  because  of  strain, 
the  confining  potential  for  each  carrier  type  is  shifted. 
Since  the  strain  varies  with  position,  the  confining 
potentials  will  also  vary  with  position.  The  material 
parameters  used  in  our  computations  are  given  in 
Table  2.  Under  these  considerations,  the  conduction 
band  can  be  given  as 

£c(s)  =  £c0  +  <52<c(£), 


(18) 


V-G.  Stoleru  et  at.  /  Physica  E  15  (2002)  131-152 


143 


vo 

c 

CO 


r — '  ’  " '  n 

Eu  : 

h 

- TPQD3 

. TPQD2 

----  TPQD1 

— j — , — , — , — , — i — j — i — i — • — i — ' — r— 

-15  -10  -5  0  5  10  15 


(d) 


z  (nm) 


Fig.  9.  Strain  tensor  components  e^  and  ezz  along  the  z-axis  for  t  =  0.25  (TPQD1),  t  =  0.50  (TPQD2),  and  t  =  0.75  (TPQD3):  (a)  Sxx, 
(b)  Ezz,  (c)  £b,  and  (d)  Eh. 


where  E®  is  the  offset  of  the  unstrained  conduction 
band,  which  is 

£c=£v°,av  +  y  +Ee,  (19) 

and  SEC  is  the  strain-induced  shift  of  the  conduction 
band  which  is  expressed  as 

SEc(e)  =  Gc(8xx  “I"  £yy  4"  &zz\  (20) 

In  Eq.  (19),  Jo  is  the  spin-orbit  splitting,  E%  is  the  un¬ 
strained  band  gap,  and  E^  av  is  the  unstrained  average 
valence  band-edge.  The  parameter  ac  in  Eq.  (20)  is 
the  deformation  potential  for  the  conduction  band. 

The  effect  of  strain  on  the  valence  band  depends 
largely  on  the  symmetry  of  the  strain.  The  heavy-  and 
light-hole  energy  bands,  is**  and  is*,  couple  to  the 
individual  strain  components  via  the  relations  [31]: 

Ef  =  <av  +  ^  +  <5£v,h  -  l-  3EVth,  (21 ) 


and 

£?  =<av  -  y  +  <5£v,  h  +  \  SEVt  b 

+  X-  \j A\  +  Ao5Ev^  +  \  (^v,b )2,  (22) 

where  5EV ?h  =  av£h>  and  <5i?v,b  —  be b-  The  deformation 
potentials  av  and  b  are  given  in  Table  2. 

The  dominant  effect  of  the  strain  is  that  the  dot  expe¬ 
riences  a  large  increase  in  its  band  gap  due  to  the  con¬ 
siderable  hydrostatic  pressure.  The  conduction  band 
for  structure  PQD3  in  Table  1  has  a  potential  well  that 
is  0.4  eV  deep  at  the  base  of  the  dot,  tapering  off  to  a 
depth  of  about  0.27  eV  at  the  apex.  The  valence  band 
has  a  more  complicated  structure.  If  we  could  some¬ 
how  turn  off  the  strain,  the  holes  would  be  confined 
to  the  InAs  by  a  well  that  is  only  85  meV  deep.  How¬ 
ever,  the  presence  of  strain  alters  this  considerably, 
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Fig.  10.  Strain  tensor  components  in  the  x-z  plane  for  structure  TPQD4:  (a)  e** ,  (b)  eyy ,  (c)  eZZy  and  (d)  eb  (strain  expressed  in  %). 


Table  2 

Material  parameters  used  in  calculations 


Parameter 

GaAs 

InAs 

InxGai  _jfAs 

a  (A) 

5.6503 

6.0553 

(5.6503  +  0.405.x)  [32,39,40] 

Et  (eV) 

1.518 

0.413 

(1.518  -  1.580x  +  0.475X2)  [30,33,39,40]  [T  =  6.4  K] 

Et  (eV) 

1.424 

0.324 

[0.324  +  0.7(1  -  x)  +  0.4(1  -  x)2]  [39,40]  [T  =  300  K] 

Ao  (eV) 

0.340 

0.380 

(0.340  -  0.093x  +  0.133x2)  [30,38,39] 

C„  (N/m2) 

11.88 

8.33 

(11.88  -  3.55x)  [32,33,39,40] 

Ci2  (N/m2) 

5.38 

4.53 

(5.38  -  0.85x)  [32,33,39,40] 

ac  (eV) 

-8.013 

-5.08 

(-8.013 +  2.933x)  [30,33,39] 

flv  (eV) 

1.16 

1.00 

[30,33,39] 

b(eV) 

-1.7 

-1.8 

[30,33,39] 

E°v, av  (eV) 

-6.92 

-6.747 

(-6.92  +  0.23 lx  -  0.058x2)  [30,33,39] 

ml 

0.0667wo 

0.02226wo 

(0.0667  -  0.04 19x  -  0.00254x2)  [33,39] 

EP  (eV) 

25.7 

22.2 

[(1.238  -  0.2095x)(l  -  m,/m,)(3Eg(Eg  +  A0)/3Eg  +  2J0)]  [33,38] 

and  it,  in  fact,  makes  the  dominant  contribution  to  the 
hole  confinement  potential.  A  remarkable  feature  of 
the  valence  band  is  that  it  is  peaked  at  the  apex  of  the 


dot  and  near  the  base.  This  is  clearly  evident  in  the 
energy  band  diagram  presented  along  the  z-axis  for 
structure  PQD3  in  Fig.  1 1 . 
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Fig.  11.  Electron  (Fe),  heavy-hole  (Fhh),  and  light-hole  (Fjh) 
potential  profiles  for  structure  PQD3,  plotted  along  the  z-axis. 


The  confining  potentials,  inclusive  of  the  effects 
of  strain,  are  piecewise  continuous  functions  of  posi¬ 
tion.  These  potentials  are  shown  in  Fig.  1 1  for  both 
electrons  and  holes.  The  split-off  valence  band  is  suffi¬ 
ciently  far  off  in  energy  from  the  heavy-  and  light-hole 
band-edges,  so  it  plays  no  role  in  the  calculations. 
Note  that  the  heavy-  and  light-hole  confining  potential 
wells  are  shown  inverted.  The  anisotropic  (biaxial) 
components  of  the  strain  in  the  dot  and  barrier  reduce 
the  symmetry  of  the  conventional  cubic  unit  cells, 
lifting  the  heavy-  and  light-hole  degeneracy.  The 
compressive  strain  in  the  barrier  shifts  the  GaAs 
conduction  band-edge  slightly  above  the  unstrained 
level  of  1.52  eV.  Note  that  in  Fig.  11,  the  light-hole 
band-edge  is  higher  in  energy  than  the  heavy-hole 
band-edge  in  the  barrier,  and  at  the  apex  of  the  pyra¬ 
mid.  The  heavy-hole  band  is  the  uppermost  band 
at  the  base  of  the  pyramid.  The  direction  and  mag¬ 
nitude  of  the  splitting  of  the  light-  and  heavy-hole 
bands — in  the  absence  of  appreciable  shear  strain 
components — is  dependent  solely  on  the  magnitude 
and  sign  of  the  biaxial  strain,  £*>•  In  those  regions  of 
the  structure  where  the  biaxial  strain  is  negative,  the 
light-hole  band  will  be  shifted  upwards  in  energy  and 
the  heavy-hole  band  downwards;  in  those  regions 
where  the  biaxial  strain  is  positive,  the  heavy-hole 
band  will  be  uppermost.  When  the  biaxial  strain  is 
zero,  the  light-  and  heavy-hole  bands  are  degenerate. 

We  show  the  electron  confining  potential  in  the 
x-z  plane  in  Fig.  12(a)  for  structure  PQD3.  Here,  the 
zero  of  energy  is  fixed  at  the  GaAs  level.  The  potential 
well  for  electrons  has  a  depth  of  about  400  meV  over 
much  of  the  pyramid,  deepening  to  about  450  meV  at 
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Fig.  12.  Confining  potential  experienced  by  charge  carriers  in  the 
x-z  plane  for  structure  PQD3  (negative  values  represent  an  attrac¬ 
tive  potential);  (a)  electrons,  (b)  heavy-holes,  and  (c)  light-holes 
(potential  expressed  in  eV). 
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the  base.  In  the  wetting  layer,  the  potential  is  taken 
to  be  identical  to  the  potential  near  to  the  base  of  the 
pyramid.  The  potential  in  the  barrier  is  close  to  zero 
since  the  material  here  is  GaAs.  Fig.  12(b)  shows  the 
contour  plot  of  the  confining  potential  experienced 
by  heavy-hole  carriers.  In  the  regions  attractive  to 
heavy-holes,  the  potential  is  negative.  The  potential 
inside  the  well  has  a  positive  gradient  from  the  base 
towards  the  apex.  Over  this  distance,  the  potential 
changes  by  several  hundred  meV.  Because  of  this  and 
the  large  effective  mass  for  heavy-holes,  one  would  ex¬ 
pect  charge  localization  near  the  base  of  the  pyramid. 
The  contour  plot  for  the  light-hole  confining  potential 
is  shown  in  Fig.  12(c).  This  plot  shows  a  slowly  vary¬ 
ing  attractive  potential  in  the  barrier  region  above  and 
below  the  pyramid.  It  reaches  a  minimum  at  the  apex 
of  the  pyramid. 

Several  methods  have  been  developed  to  calculate 
the  electronic  structure  of  the  square-based  pyramidal 
InAs  dots  [4,6,11].  Among  these  is  the  pseudopoten¬ 
tial  plane-wave  approach  used  by  Williamson  et  al. 
[32].  To  calculate  the  energy  levels  and  electron  (or 
hole)  wave  functions  we  use  an  eight-band  effective 
mass  approach.  The  strain  effect  is  included  via  defor¬ 
mation  potential  theory  [33].  The  bound  energy  levels 
are  computed  as  functions  of  quantum-dot  size.  The 
bound  states  of  the  dot  are  found  by  numerically  solv¬ 
ing  the  Schrodinger  equation,  which  in  the  effective 
mass  approximation,  can  be  written  as 

-Tv(^))v,'-W+,'W5,“W 

=EV„(r).  (23) 

In  this  expression,  m*{r)  must  be  replaced  by 
wfnGaAsO)  inside  the  dot,  and  by  m*GllAs(r)  in  the 
matrix  material;  V(r)  is  the  three-dimensional  con¬ 
fining  potential.  In  the  framework  of  the  eight-band 
model,  the  wave  function  can  be  expanded  as  a  linear 
combination  of  the  basis  functions,  thus 

8 

'Pn(r)  =  YsWujir)  =  (24) 

j=\  3,JZ 

where  uj(r)  has  the  periodicity  of  the  crystal  lattice 
and  j  is  the  band  index.  At  the  band-edges,  these 
functions  are  characterized  by  symmetry  arguments  as 
eigenstates  \J,J2)  of  the  Bloch  angular  momentum  J. 


The  states  \J,JZ)  are  the  band-edge  functions  of  the 
Bloch  state  space  [34],  These  wave  functions  are,  for 
the  ^-like  T\  conduction  band, 

I,I}=|5)|T)  and 

(25) 

while  for  the  y>-like  Tg,  T 7  valence  bands,  they  can 
be  written  as  [34] 

2'  2X  =  (V2)  <|jr>  +  T>’ 

5'0,=(^)<W+'1>,>)U> 

-  (If)  ™ 

-  (vf) |z)l 

i’“5)v=(yi)(W",W)U)' 

?i),=(ys)<w+'w)ll> 

+  (If)  »i  T>- 

+  (\ff)  Wl  i>-  m 

We  want  to  point  out  that  it  is  generally  known  that 
the  k  p  method  has  some  problems  associated  with  it 
when  applied  to  calculations  involving  nanostructures 
[33].  These  difficulties  include  (i)  the  fixed  number  of 
Bloch  functions  (eight  in  the  formalism  used  here), 
(ii)  the  restriction  of  the  validity  of  the  method  to  the 
Brillouin  zone  center,  (iii)  the  use  of  the  same  Bloch 
functions,  regardless  of  material  and  strain  variations, 
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and  (iv)  the  difficulty  of  choosing  appropriate  bound¬ 
ary  conditions  with  associated  matching  criteria  for 
the  envelope  functions  across  the  heterointerfaces.  In 
our  calculations,  we  take  into  account  the  variation 
of  mass  parameters  from  their  bulk  values  due  to  the 
strain-induced  band  deformations.  The  details  of  the 
guidelines  used  have  been  discussed  by  Burt  [35]. 

The  eigenvalues  and  eigenfunctions  of  Eq.  (23)  are 
obtained  by  invoking  periodic  boundary  conditions, 
expanding  'Fn(r)  in  terms  of  normalized  plane-wave 
states,  and  diagonalizing  the  resulting  matrix.  This 
approach,  which  has  also  been  used  by  Cusack  et  al. 
[6],  does  not  require  matching  of  the  wave  functions 
across  the  boundary  between  the  dot  and  the  matrix 
materials.  This  makes  the  method  applicable  to  arbi¬ 
trary  (position  dependent)  confining  potentials.  The 
only  requirement  on  the  boundary  conditions  is  that 
the  states  do  not  significantly  overlap  for  dots  adjacent 
to  each  other. 

In  general,  the  number  of  confined  states  in  a  quan¬ 
tum  dot  depends  on  the  size  of  the  dot  and  on  the  thick¬ 
ness  of  the  wetting  layer.  For  the  conduction  band, 
there  are  usually  only  a  few  bound  states.  The  energy 
spacing  between  the  ground  state  and  the  first  excited 
state  in  the  conduction  band  typically  ranges  from 
about  60  to  95  meV  for  the  dot-sizes  considered  in 
our  work  (see  Fig.  13).  We  find  that  the  valence-band 
states  are  more  tightly  confined  because  of  the  large 
hole  effective  mass.  The  energy  spacing  here  ranges 
from  a  few  meV  to  30  meV. 

In  our  calculations,  we  have  neglected  the  Coulomb 
interaction  energy.  This  energy  largely  depends  on 
the  value  of  the  dielectric  constant.  As  a  result,  dots 
of  the  same  size  can  belong  to  different  confinement 
regimes  in  materials  with  different  dielectric  constants. 
In  III-V  compounds,  with  a  typical  relative  dielectric 
constant  of  ~13,  the  bulk  exciton  radius  is  >  10  nm, 
causing  a  structural  quantum  dot  of  similar  dimension 
and  sufficiently  deep  potential  to  be  in  the  strong  con¬ 
finement  regime  [22].  Because  of  this,  additional  bind¬ 
ing  energy  from  the  Coulomb  effect  is  negligible  since 
the  dots  are  already  in  the  strong  confinement  regime. 
In  any  case,  the  Coulomb  interaction  energy  is  small 
compared  to  the  separation  of  the  ground  state  from 
the  first  excited  state  for  both  electrons  and  holes. 

In  Fig.  13,  we  show  the  calculated  electron  and  hole 
energy  level  dependence  on  pyramid  base  length.  The 
electron  (hole)  levels  are  plotted  relative  to  the 


Fig.  13.  (a)  Electron  and  (b)  hole  quantum-dot  energy  levels 
(ground  state  and  first  two  excited  states),  displayed  as  a  function 
of  dot  base  size. 

unstrained  GaAs  conduction  (valence)  band-edge.  For 
dot  base  dimensions  smaller  than  nm,  no  bound 
electron  states  are  predicted.  Beyond  6nm,  a  few 
states  are  predicted  in  the  conduction  band.  And  as 
stated  earlier,  many  states  are  predicted  in  the  valence 
band.  This  is  due  to  the  large  effective  mass  associ¬ 
ated  with  the  holes  and  to  the  nature  of  the  light-hole 
confining  potential  whose  smoothly  varying  form 
leads  to  a  quasi-continuum  of  tenuously  bound  states. 

We  have  also  determined  the  envelope  functions  for 
the  first  few  confined  states  in  the  dot.  Fig.  14(a),  for 
example,  shows  the  modulus-squared  envelope  func¬ 
tion  for  the  Eq o  state  for  structure  PQD3;  this  function 
is  plotted  in  the  y-z  plane,  cutting  through  the  pyramid 
and  the  wetting  layer.  The  relatively  isotropic  char¬ 
acter  of  the  confining  potential  for  electrons,  coupled 
with  the  small  effective  mass,  results  in  a  state  that 
permeates  throughout  the  dot  and  penetrates  into  the 
sides  of  the  pyramid.  Fig.  14(b)  is  the  hole  envelope 
function  in  the  y-z  plane  for  the  state  E\&.  Unlike  the 
ground  state  in  the  conduction  band,  the  hole  ground 


148 


V.-G.  Stoleru  et  aUPhysica  E  15  (2002)  131-152 


75: 

5  \ 

2.5 : 

-2.5: 

JjfjfHI  'iBBk 

-5: 

-75; 

(c)  -10  -7.5  -5  -2.5  0  2.5  5  75  10 


F|2  z(nm) 


Above 

0.99  < 

Hi 

0.90 

0.99 

B 

0.81 

0.90 

H 

0.77 

0.81 

p 

0.63 

0.77 

||| 

0.54 

0.63 

9 

0.45 

0.54 

0.36 

0.45 

0.27 

0.36 

0.18 

0.27 

"i 

0.09 

0.18 

\ 

Below 

0.09 

(am) 


Fig.  14.  The  squared  absolute  magnitude  of  the  envelope  function  for  structure  PQD3  across  the  y-z  plane  for  energy  levels:  (a)  £eo,  (b) 
£ho,  (c)  EhU  (d)  Ev 


state  is  confined  to  the  base  of  the  dot.  In  Fig.  14(c) 
and  (d),  we  show  the  envelope  functions  for  the  hole 
excited  states  2?hi  and 

The  calculated  transition  energies  for  the  ground 
state,  Eqo  — *  jEhO>  and  for  the  first  £ei  —►-Em*  and  sec¬ 
ond  Et 2  — >  E\a  excited  states  for  pyramids  of  various 
base  dimensions  are  shown  in  Fig.  15(a).  From  Fig. 
15(a),  the  calculated  fundamental  transition  energy  for 
a  pyramid  with  a  base  length  of  16  nm  is  1.10  eV. 
This  value  is  in  good  agreement  with  the  experimen¬ 


tally  determined  transition  peak  at  1.098  eV  shown 
in  Fig.  15(b).  The  same  good  agreement  is  found 
between  the  calculated  and  measured  peak  lumines¬ 
cence  energy  values  for  a  20-nm-base  pyramidal  quan¬ 
tum  dot;  the  measured  experimental  data  is  shown  in 
Fig.  15(c). 

The  energy  splitting  between  the  ground  state  and 
first  excited  hole  state  of  30  meV  in  Fig.  15(b)  is  in 
good  agreement  with  a  recent  experimental  study  of 
the  sub-level  structure  which  measured  a  difference 
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Fig.  15.  (a)  Calculated  transition  energies  for  pyramidal  QDs: 
Ec0-Eh0,  £ei-iThi>  and  £e2-£h2,  as  functions  of  QD  base  width; 
(b)  Fundamental  (£eo-£ho)  transition  energy  measured  by  PL  for 
a  dot  of  base  B—\6  nm;  and  (c)  for  a  dot  of  base  width  5=20  nm. 


of  approximately  27  meV  [36].  The  Eeo  — ►  £ho 
transition  is  the  dominant  excitation  in  all  of  the 
structures  studied. 

We  now  consider  the  computation  of  the  electronic 
energy  levels  in  truncated  pyramidal  quantum  dots. 
The  basic  approach  is  similar  to  what  has  already 
been  discussed  in  this  paper.  For  truncation  factors 
ranging  from  0.25  to  0.75,  we  show  the  confining 
potentials  for  electrons  and  holes  in  Fig.  16.  The 
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Fig.  16.  Confining  potentials  for  electrons  and  heavy-holes,  along 
z-axis,  for  structures  TPQD(l-3)  with  truncation  factors:  /  =  0.25 
(TPQD1),  t  =  0.50  (TPQD2),  and  t  =  0.75  (TPQD3). 


bound  states  are  determined  by  numerically  solv¬ 
ing  the  Schrodinger  equation  within  the  context  of 
approximations  similar  to  those  used  in  the  ideal 
pyramid  case.  The  calculated  energy  levels  for  the 
ground  state  and  two  excited  states  in  the  conduction, 
as  well  as  the  valence  bands  are  shown  in  Fig.  17  for 
a  range  of  truncation  factors.  As  a  specific  example, 
we  calculated  the  energy  levels  for  structure  TPQD3, 
whose  dimensions  are  given  in  Table  1;  the  TEM 
micrograph  of  this  structure  was  shown  in  Fig.  8. 
The  first  three  confined  electron  energy  levels,  mea¬ 
sured  with  respect  to  the  unstrained  GaAs  conduction 
band-edge,  are  Et o  =  1.404  eV,  EQ\  =  1.426  eV,  and 
Eq2  =  1.428  eV.  In  the  valence  band,  the  correspond¬ 
ing  heavy-hole  levels — again  measured  with  respect 
to  the  unstrained  GaAs  valence  band-edge — are 
£ho  =  328.79  meV,  E^\  =  343.88  meV,  and  Ey&  — 
345.74  meV.  The  transition  energies,  for  the  allowed 
transitions,  are  shown  in  Fig.  18(a);  the  fundamental 
transition,  Et o  — >  £ho,  is  at  1.076  eV.  This  energy 
is  in  good  agreement  with  the  experimentally  deter¬ 
mined  peak  of  1.08  eV  for  the  photoluminescence 
emission  spectrum  shown  in  Fig.  18(b)  for  structure 
TPQD3.  The  photoluminescence  measurements  were 
carried  out  using  an  Ar+  ion  laser  (2  =  488  nm)  as  an 
excitation  source.  The  emitted  radiation  was  detected 
with  a  liquids -cooled  Ge  detector. 

In  devices  such  as  near-infrared  lasers  and 
mid-infrared  detectors,  the  preferred  medium  in  the 
active  region  is  often  an  (In,Ga)As/GaAs,  rather  than 
an  InAs/GaAs  quantum-dot  superlattice.  In  this  case, 
it  becomes  necessary  to  perform  the  energy  level 
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Fig.  17.  Ground  state  and  first  two  excited  states  for  electrons 
and  heavy-holes  for  structures  TPQD(l-3)  with  truncation  factors: 
t  =  0.25-0.75. 


calculations  for  (In,Ga)As  dots.  We  have  per¬ 
formed  such  calculations  for  truncated  pyramidal 
Ino.23Gao.77 As  dots.  The  parameters  used  in  the  calcu¬ 
lations  are  given  in  Table  2.  For  illustrative  purposes, 
we  have  used  the  truncated  structure  TPQD4  in  Table 
1 .  The  fundamental  transition  energy,  2seo  — 5 ►  £h0,  for 
such  a  structure  is  calculated  to  be  1.126  eV.  In  the 
conduction  band,  the  electron  ground  state  of  such 
a  dot  is  separated  by  about  107  meV  from  the  first 
excited  state.  These  computed  values  are  in  good 
agreement  with  the  photoluminescence  and  infrared 
absorption  data  obtained  at  300  K  [37]. 

5.  Conclusions 

The  relaxation  of  strain  and  its  residual  compo¬ 
nent  in  lattice-mismatched  epitaxy  is  responsible 
for  the  self-organization  of  quantum  dots  in  the 
Stranski-Krastanow  crystal  growth  mode.  Beginning 


Fig.  18.  Ground  state  and  first  excited  state  transition  energies 
for  structures  TPQD(l-3)  with  truncation  factors:  t  =  0.25-0.75: 
(a)  calculated,  (b)  measured. 


from  this  premise,  we  have  calculated  the  strain 
distributions  in  pyramidal  and  truncated  pyramidal 
(In,Ga)As/GaAs  quantum  dots  using  a  method  based 
on  Eshelby’s  inclusion  theory  of  continuum  elastic¬ 
ity.  It  is  found  that  the  hydrostatic  component  of  the 
strain  is  mostly  confined  within  the  dots,  while  the 
biaxial  strain  is  transferred  from  the  dot  to  the  barrier 
material.  By  taking  into  account  the  influence  of  the 
strain  on  the  band  gap  of  the  dots,  we  have  solved 
the  three-dimensional,  effective  mass,  single-particle 
Schrodinger  equation  for  the  electronic  energy  levels 
in  the  dot.  We  find  that  strain  plays  a  major  role  in 
the  energy  structure  of  the  quantum  dots.  For  the 
pyramidal  geometry  considered  here,  the  electronic 
energy  levels  are  also  a  sensitive  function  of  the  base 
length. 

The  results  of  our  calculations  are  in  good  agree¬ 
ment  with  those  reported  in  the  literature,  even 
though  our  method  is  considerably  simpler  that  those 
used  by  others  [4,20].  For  the  fundamental  transition 
of  a  ground  state  electron  in  the  conduction  band 
recombining  with  a  hole  in  its  ground  state  in  the 
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valence  band,  we  find  that  our  experimental  results 
also  agree  with  calculations. 

In  summary,  we  have  developed  a  simple  method 
for  rapidly  calculating  strain  in  embedded  quantum 
dots.  The  method  allows  one  to  determine  the  inter¬ 
dependence  of  strain  on  shape,  and  by  extension,  the 
effects  of  strain  on  confined  energy  levels  in  quantum 
dots. 
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Appendix  A 

This  section  gives  the  more  general  mathematical 
expressions  for  the  stress  distributions  inside  pyrami¬ 
dal  InAs  quantum  dots  of  different  degrees  of  trunca¬ 
tion.  The  expressions  are  obtained  by  integrating 
Eq.  (6).  They  are  given  as  functions  of  the  space  co¬ 
ordinate  z,  in  the  growth  direction.  The  strain  compo¬ 
nents  follow  immediately  from  Hooke’s  law,  as  given 
in  Eq.  (15).  The  following  definitions  are  used  in  the 
expressions:  the  pyramid  base  width  is  B  =  2a;  the 
height,  in  the  absence  of  truncation,  is  H ;  the  param¬ 
eter  t  is  the  degree  of  truncation;  and  h  is  the  height 
of  a  truncated  pyramid.  The  stress  component  oxx,  for 
example,  is  written  as  a  summation  over  three  terms: 

_  h  Eb o 

Gxx  a 2  +  h 2  47t(1  —  v) 

X  [°xt(l)  +  &xx(2)  +  °xt(3)]>  (A- 1 ) 

where  the  terms  o^(2)>  an<l  a'xx(i)  are  giyen  as 
,  4a2  f,  \-2a\\-t)-h(z-ht) 

"  7WW  l  [  vPtf 
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Similarly,  for  ayy  we  derive  the  following 
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The  other  stress  component,  <r2Z,  follows  from: 
2Eeq 
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Abstract 

This  paper  reports  on  interband  emission  and  intraband  absorption  processes  in  doped  InGaAs/GaAs  quantum-dot  nano¬ 
structures  as  the  measurement  temperature  is  varied.  An  analysis  of  the  data  yields  useful  insight  into  the  carrier  redistribution 
mechanisms  that  occur  as  a  result  of  dot-size  variation.  For  the  samples  studied,  the  size  distribution  is  bimodal;  this  is  reflected 
in  the  results  obtained  from  the  interband  emission,  as  well  as  in  the  intraband  absorption  data.  The  understanding  gained 
from  studies  such  as  these  is  important  in  the  design  of  quantum-dot  optoelectronic  devices  intended  for  room-temperature 
operation.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 

PACS:  73.20.Dx;  78.55.Cr;  78.66.Fd 

Keywords:  Quantum  dots;  Photoluminescence;  Absorption;  Photodetectors;  InGaAs/GaAs 


Zero-dimensional  quantum-dot  nanostructures  have 
recently  become  of  considerable  interest  because  of 
their  potential  application  in  electronic  and  optoelec¬ 
tronic  devices.  Some  of  the  important  characteristics 
of  the  nanostructures  include:  atom-like  density  of 
states,  large  exciton  binding  energies,  and  tunability 
of  the  emission  and  absorption  wavelengths.  These 
unique  characteristics  have  stimulated  a  number  of 
efforts  aimed  at  developing  a  new  generation  of 
light-emitters  and  detectors  [1,2]. 


*  Corresponding  author.  Fax:  +1-434-924-8818. 

E-mail  addresses:  dp4b@virginia.edu  (D.  Pal), 
towe@cmu.edu  (E.  Towe). 

1  On  leave  from  the  St.  Petersburg  State  Technical  University, 
St.  Petersburg,  195251,  Russia. 

2  Present  address:  Carnegie  Mellon  University,  Department  of 
Electrical  and  Computer  Engineering,  5000  Forbes  Avenue,  Pitts¬ 
burgh,  PA  15213-3891,  USA. 


In  general,  the  ideal  electronic  or  optoelectronic 
device  should  operate  at  room  temperature.  To  en¬ 
able  this,  it  is  desirable  to  have  a  clear  understand¬ 
ing  of  the  effects  of  temperature  on  the  electronic 
and  optical  properties  of  the  material  out  of  which 
the  device  is  fabricated.  There  have  been  some 
reports  on  the  optical  properties  of  quantum-dot 
nanostructures;  specifically,  there  have  been  reports 
on  the  temperature-dependence  of  interband  emis¬ 
sion  from  InAs/GaAs  [3],  (In,Ga)As/GaAs  [4]  and 
(Al,In)As/(Al,Ga)As  [5]  quantum-dot  structures.  The 
emission  from  these  structures  has  been  observed  to 
decrease  with  increasing  temperature.  This  is  prob¬ 
ably  due  to  the  ejection  of  carriers  from  the  dots 
by  thermal  excitation  into  the  lowest-laying  states 
in  the  wetting  layer  or  into  similar  states  in  the  ma¬ 
trix  material.  When  this  occurs,  the  usual  radiative 
recombination  processes  no  longer  take  place  within 
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the  dot — leading  to  a  decrease  in  emission.  A  second 
possible  pathway  is  a  non-radiative  recombination 
process;  yet  another  pathway  for  quenching  is  the 
thermal  redistribution  of  carriers  within  dots  of  dif¬ 
ferent  sizes  [6].  The  temperature-dependent  carrier 
redistribution  within  dots  of  different  sizes  is  gener¬ 
ally  thought  to  be  responsible  for  the  unusual  line 
width  variation.  Furthermore,  it  is  connected  with  the 
rapid  red-shift  of  the  photoluminescence  emission 
peak.  This  red-shift  is  more  pronounced  than  the  usual 
temperature-dependent  band  gap  change.  We  want 
to  point  out  that  the  thermal  redistribution  of  carri¬ 
ers  in  doped  quantum-dot  structures— such  as  those 
used  in  quantum-dot  infrared  photodetectors — is  not 
yet  fully  understood.  The  redistribution  of  carriers  at 
high  temperatures  can  modify  the  intraband  absorp¬ 
tion  on  which  these  detectors  depend.  This,  in  turn, 
can  affect  the  performance  of  the  devices  at  elevated 
temperatures. 

Normal-incidence  intraband  absorption  has  been 
reported  in  n-type  doped  InGaAs/GaAs  [7]  and 
InAs/GaAs  [8]  quantum  dots.  Phillips  et  al.  have 
reported  intraband  absorption  in  a  five-period 
InAs/GaAs  quantum-dot  structure  at  300  K  [9].  Other 
researchers  have  reported  normal-incidence  intraband 
absorption  processes  in  InAs/InALAs  quantum-dot 
structures  at  77  K  [10].  However,  there  has  not  been 
any  reports  on  the  effect  of  carrier  redistribution  on 
the  absorption  strength. 

In  this  paper,  we  present  and  discuss  temperature- 
dependent  results  for  interband,  emission,  as  well 
as  normal-incidence  intraband  absorption  in  doped, 
five-period  InGaAs/GaAs  quantum-dot  nanostruc¬ 
tures.  The  dot-sizes  in  these  structures  fall  roughly 
into  two  groups:  one  group  has  a  large  dot-size,  while 
another  has  a  small  dot-size.  We  find  that  the  relative 
populations  of  the  states  in  the  two  dot  groups  follow 
Boltzmann  statistics.  Finally,  we  measure  an  intra¬ 
band  absorbance  of  about  7.0%  in  our  quantum-dot 
structures  at  77  K. 

Our  structures  were  grown  on  (0  0  1)  GaAs  sub¬ 
strates  in  a  solid-source  molecular-beam  epitaxy 
system.  A  typical  structure  consists  of  five  periods 
of  InGaAs/GaAs  quantum-dot  superlattice  layers. 
The  nominal  Indium  composition  based  on  growth 
parameters  and  in  situ  RHEED  observations,  was 
about  30%.  However  the  X-ray  analysis  of  the  sample 
showed  that  the  average  Indium  composition  in  the 


Fig.  1.  AFM  image  (1  Jim  x  1  ^im)  of  InGaAs  quantum  dots. 


dots  is  about  23%.  The  dots  were  delta-doped  with  Si 
to  yield  about  two  electrons  per  dot.  The  thickness  of 
the  GaAs  barrier  layer  in  the  superlattice  was  about 
50  nm.  For  atomic  force  microscope  (AFM)  studies, 
several  uncapped  samples  were  grown.  These  samples 
were  grown  under  conditions  identical  to  those  used 
during  the  growth  of  the  experimental  structures. 

Fig.  1  is  a  1  pm  x  1  pm  AFM  micrograph  of  one 
of  the  uncapped  InGaAs  samples.  The  dot  density 
is  ~  6.34-6.8  x  1010  cm-2.  Notice  that  the  dots  are 
closely  packed  and  are  not  uniform  in  size.  The  dot 
heights  vary  from  about  5-7.5  nm.  We  have  plot¬ 
ted  the  statistical  distribution  of  the  lateral  dot-size  in 
the  histogram  shown  in  Fig.  2.  This  histogram  shows 
a  bimodal  dot-size  distribution.  The  average  lateral 
extent  of  the  majority  of  the  dots  is  either  ~  28  or 
~  34  nm,  as  shown  in  the  histogram.  Note  that  there  are 
apparently  more  small-size  dots  than  large-size  dots. 

Photoluminescence  measurements  were  carried 
out  on  the  samples  in  the  temperature  range  from 
77  to  300  K  using  an  Argon  ion  laser  (488  nm)  as 
the  excitation  source.  The  emitted  radiation  was  de¬ 
tected  with  a  liquid-nitrogen-cooled  Ge  photodetector. 
Normal-incidence  intraband  absorption  measurements 
were  performed  with  a  BOMEM  Fourier  Transform 
Infrared  spectrometer  in  the  temperature  range  from 
77  to  300  K. 
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Fig.  2.  Histogram  of  the  lateral  size  distribution  of  InGaAs  dots. 


Fig.  3  shows  typical  photoluminescence  emission 
spectra  from  one  of  our  samples.  At  the  low  tem¬ 
perature  of  77  K,  the  emission  is  characterized  by  a 
single,  narrow  peak  at  1.19  eV,  with  a  spectral  spread 
of  about  30.8  meV.  As  the  temperature  is  increased, 
the  spectrum  broadens,  and  at  about  300  K,  in  addi¬ 
tion  to  the  main  peak  at  1.12  eV,  there  is  a  distinct 
identifiable  shoulder  on  the  high-energy  side  of  the 
spectrum.  A  deconvolution  of  this  spectrum  reveals 
the  two  spectra  labeled  (a)  and  (b),  with  peak  ener¬ 
gies  at  about  1.12  and  1.16  eV,  respectively.  The  two 
peaks  are  separated  by  about  40  meV;  this  separation 
is  smaller  that  the  expected  energy  difference  of  80- 
100  meV  between  the  ground  state  and  the  first  ex¬ 
cited  state  in  the  conduction  band  of  our  typical  dots. 
From  intraband  absorption  measurements — which  we 
will  discuss  later — the  separation  between  the  ground 
state  and  the  first  excited  state  in  the  conduction  band 
is  about  100  meV.  This  information  suggests  that  the 
emission  peak  at  1.16  eV  is  probably  not  related  to 
an  interband  transition  from  the  first  excited  state 
in  the  conduction  band  of  a  certain  group  of  dots 
to  the  relevant  energy  level  in  the  valence  band.  A 
possible  explanation  is  that  the  emission  is  associ¬ 
ated  with  an  interband  transition  from  the  ground 


Fig.  3.  The  low  (77  K)  and  room  temperature  (300  K)  photolu¬ 
minescence  spectra  of  a  five-period  (In,Ga)As/GaAs  quantum-dot 
nanostructure. 

state  of  a  different  set  of  dots  with  a  different  average 
size. 

The  appearance  of  a  suggestive  shoulder  (and  hence 
a  second  peak)  in  the  emission  spectrum  may  be 
explained  if  we  consider  the  bimodal  dot-size 
distribution  observed  in  the  AFM  images  of  the  sam¬ 
ples.  The  bimodal  size  distribution  means  that  the 
dot  population  falls  into  two  clusters,  with  the  av¬ 
erage  size  of  one  set  of  dots  clearly  identifiable  as 
“small”,  and  another  set  with  an  average  size  that  is 
“large”.  During  the  optical  excitation  process,  carri¬ 
ers,  particularly  holes  in  doped  dots,  tunnel  out  of  the 
“small”  dots  into  the  “large”  ones  so  that  they  can 
occupy  the  lowest  possible  energy  states  available 
in  the  ensemble.  At  77  K,  the  luminescence  from 
the  recombination  radiation  is  therefore  primarily 
from  the  dots  that  have  the  lowest  possible  ground 
states — that  is:  the  large-size  dots.  This  is  the  lu¬ 
minescence  whose  peak  is  observed  at  1.19  eV.  As 
the  temperature  is  increased,  thermal  excitation  be¬ 
gins  to  eject  carriers  out  of  the  low  energy  states, 
leading  to  a  redistribution  such  that  dots  of  all  sizes 
are  occupied.  For  the  bimodal  size  distribution  ob¬ 
served  in  the  atomic  force  microscope  images  of 
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Fig.  4.  The  relative  energy  shift,  A£qd,  as  a  function  of  tem¬ 
perature  for  bands  (a)  and  (b),  measured  with  reference  to 
an  emission  peak  at  107  K.  The  solid  line  shows  the  relative 
temperature- dependence  of  the  band  gap  energy,  A Eg,  for  InGaAs. 
The  Stokes  shift  as  a  function  of  temperature  is  also  shown  (on 
the  right-hand  side  abscissa). 

our  samples,  the  expected  room-temperature  lumi¬ 
nescence  is  shown  in  Fig.  3;  this  spectrum  can  be 
resolved  into  the  spectra  (a)  and  (b).  In  this  figure, 
spectrum  (a)  is  attributed  to  transitions  associated 
with  the  large-size  dots;  this  is  because  the  energy 
levels  are  located  deeper  at  the  bottom  of  the  potential 
wells  of  the  dots.  The  peak  at  1.16  eV,  in  spectrum 
(b),  is  a  result  of  transitions  within  the  small  dots. 

The  thermal  redistribution  of  carriers  can  be 
followed  systematically  by  observing  a  relative 
shift,  AEqo — from  a  fixed  reference  point — for 
each  of  the  two  emission  peaks  as  the  temperature 
is  changed.  We  have  chosen  the  single  emission 
peak,  ftvio7K>  observed  at  107  K  as  the  fixed  ref¬ 
erence  point.  The  relative  shift  can  then  be  given 
as  A£qd  =  hv(T)  -  Av(107K),  where  hv(T)  is 
the  temperature-dependent  peak.  Fig.  4  illustrates 
the  relative  temperature-dependent  shift  for  the  two 
spectra  (a)  and  (b)  of  Fig.  3.  For  comparison,  we 
show  (with  the  solid  line)  the  relative  shift  of  the 
band  gap  energy  of  Ino.23Gao.77As  as  calculated  from 
AEg  =  Eg(T)  -  £g(107  K)  according  to  Varshini’s 
law  [11]. 

In  general,  with  no  carrier  redistribution  processes 
involved,  one  would  expect  the  photoluminescence 


Fig.  5.  The  Arrhenius  plot  of  the  ratio  of  the  integrated  photo- 
luminescene  intensity  of  band  (a)  to  band  (b)  as  a  function  of 
temperature. 


emission  peak  related  to  quantum-dot  emission  to 
follow  a  trend  similar  to  that  of  the  band  gap  with 
increasing  temperature.  This,  however,  is  not  the  case. 
The  relative  rate  of  shift,  0A7sqd/9T,  is  small.  If  it 
were  larger,  it  would  indicate  carrier  redistribution 
from  high  to  low  energy  states.  This  has  been  ob¬ 
served  by  Polimeni  et  al.  in  InAs  quantum  dots  [12]. 
In  our  case,  since  the  dots  are  very  tightly  packed, 
most  of  the  carriers  initially  occupy  the  lowest  energy 
states  due  to  tunneling  processes. 

The  energy  difference,  A 2sqD  —  AEg,  is  often  de¬ 
fined  as  the  Stokes  shift  (SS)  between  the  lumines¬ 
cence  and  absorption  lines.  From  Fig.  4,  we  see  that 
the  Stokes  shift  for  spectrum  (b)  is  larger  than  that  for 
spectrum  (a).  This  suggests  that  the  dot-size  distribu¬ 
tion  corresponding  to  emission  band  (b)  is  larger  than 
that  for  band  (a). 

Another  way  to  gain  further  insight  into  the  carrier 
redistribution  processes  at  high  temperatures  is  to  plot 
the  ratio  of  the  integrated  photoluminescence  intensity 
of  spectrum  (a)  to  spectrum  (b)  as  a  function  of  tem¬ 
perature.  This  ratio  is  shown  in  Fig.  5.  The  slope  of  the 
line  in  this  figure  yields  an  activation  energy  of  about 
38  meV.  The  energy  difference  (of  40  meV)  between 
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Fig.  6.  Normal-incidence  intraband  absorbance  spectra  obtained  at 
77  and  300  K.  Spectra  (a)  and  (b)  were  obtained  by  deconvolu¬ 
tion  of  the  low-temperature  (77  K)  spectrum;  spectra  (c)  and  (d) 
were  obtained  as  a  result  of  deconvolving  the  room-temperature 
absorbance  spectrum.  A  simple  energy  band  diagram  for  the  small- 
and  large-size  dots  is  in  inset  (f)  for  the  small  dots,  and  in  inset 
(e)  for  the  large  dots. 


the  peak  of  spectrum  (a)  and  spectrum  (b)  in  Fig.  4  is 
quite  close  to  the  activation  energy;  this  lends  further 
support  to  the  notion  that  carriers  transfer  from  large- 
to-small-size  dots.  Similar  results  have  been  observed 
in  undoped  quantum-dot  structures  where  the  transfer 
of  carriers,  however,  was  from  a  dot  layer  to  the  wet¬ 
ting  layer  [13]. 

We  now  turn  to  the  discussion  of  intraband  absorp¬ 
tion  in  our  quantum-dot  nanostructures.  The  absorp¬ 
tion  spectra,  normalized  to  a  background  spectrum 
of  a  polished  semi-insulating  GaAs  substrate  without 
quantum  dots  on  it,  are  shown  in  Fig.  6.  The  high¬ 
est  absorbance  of  7%  was  recorded  for  a  sample  at 
77  K;  this  is  shown  in  spectrum  (a).  Our  absorption 
measurements  were  performed  at  normal-incidence 
on  a  sample  containing  only  five  planes  of  quantum 
dot  layers.  The  absorbance  per  quantum  dot  layer  is 
therefore  1.4%.  To  our  knowledge,  this  represents  the 
highest  absorbance  ever  reported  for  these  structures. 


We  must,  however,  comment  that  very  little  data  has 
been  reported  in  the  literature  on  absorbance.  The  few 
results  that  exist  include  data  reported  by  Weber  et 
al.  [10].  This  group  reported  low-temperature  (77  K) 
absorbance  data  for  InAs/InAlAs  quantum-dot  struc¬ 
tures.  Their  absorbance  per  dot  layer  was  0.78%.  In 
a  paper  by  Phillips  et  al.  [9],  a  room-temperature  ab¬ 
sorbance  of  about  0.48%  per  InAs/GaAs  quantum-dot 
layer  was  reported.  This  value  is  lower  than  our 
room-temperature  absorbance  of  1.02%  per  layer  of 
InGaAs  dots. 

The  low-  and  room-temperature  spectra  in 
Fig.  6  show  that  the  infrared  absorption  occurs  at 
two  energies — consistent  with  the  photolumines¬ 
cence  results  discussed  earlier.  A  deconvolution  of 
the  low-temperature  absorption  spectrum  reveals  two 
resonances:  one  at  105  meV,  as  shown  in  curve  (a), 
and  another  at  137  meV,  shown  in  spectrum  (b).  The 
room-temperature  spectrum  can  also  be  resolved  into 
two  spectra  with  peak  energies  at  107  and  135  meV. 
The  corresponding  integrated  absorbance  for  spec¬ 
trum  (a)  and  spectrum  (c)  decreases  from  3.2  x  10-3 
to  2.4  x  10-3  as  the  temperature  is  raised  from 
77  to  300  K.  On  the  other  hand,  the  integrated  ab¬ 
sorbance  increases  from  6.3  x  1CT4  for  spectrum  (b) 
to  8.2  x  10"4  for  spectrum  (d)  as  the  temperature  is 
raised  from  77  to  300  K. 

It  is  unlikely  that  the  absorption  spectra  in  Fig.  6  are 
a  result  of  transitions  from  the  ground  state  to  the  first- 
and  second-excited  states  of  the  same  set  of  dots.  This 
is  so  because  the  integrated  absorbance  of  the  relevant 
sets  of  spectra  does  not  decrease  due  to  thermal  escape 
of  carriers  from  the  dots.  The  integrated  absorbance 
decreases  only  for  the  peak  at  ~  105  meV  while  for 
the  peak  at  ~  135  meV,  it  increases  with  temperature. 

The  observed  spectra  are  more  likely  to  be  from 
intraband  transitions  in  large-  and  small-size  dots.  On 
the  basis  of  our  photoluminescence  data  and  the  intra¬ 
band  absorption  results,  one  can  derive  a  simple  pic¬ 
ture  of  the  energy  band  structure  (neglecting  Coulomb 
charging  effects)  for  the  two  sets  of  dots  under  dis¬ 
cussion.  Such  a  band  structure  is  shown  in  the  inset 
of  Fig.  6.  Here,  we  assume  that  the  energy  difference 
between  the  ground  states  in  the  valence  band  of  the 
large-  and  small-size  dots  is  approximately  2  meV.  In 
the  conduction  band,  the  energy  levels  are  as  shown. 

We  have  performed  some  theoretical  calculations 
to  determine  the  energy  levels  in  the  dots  using  an 


D.  Pal  et  al  l  Physica  E  15  (2002)  6-12 


11 


eight-band  k  •  p  model  with  the  strain  effects  included 
via  deformation  potential  theory.  The  analytical  ap¬ 
proach  developed  by  our  group  for  the  determination 
of  strain  fields  and  the  associated  energy  structure  has 
been  presented  elsewhere  [14].  The  energy  level  cal¬ 
culations  assume  that  the  large  dots  have  a  lateral 
extent  of  about  19.7  nm  and  a  height  of  3.0  nm;  the 
small  dots  have  a  height  of  3.0  nm  and  a  lateral  extent 
of  ~  16.8  nm.  The  geometry  of  the  dots  is  truncated 
pyramidal.  The  calculated  energy  levels  agree  quite 
well  with  the  experimental  results.  The  lateral  size  and 
height  of  the  dots  used  in  the  calculations  are  not  the 
same  as  those  obtained  from  the  AFM  studies.  This  is 
because:  (i)  dot-sizes  cannot  be  accurately  determined 
from  AFM  studies  due  to  finite  size  and  shape  of  the 
scanning  probe,  and  (ii)  the  lateral  size  and  height 
used  in  the  calculations  are  for  buried  dots.  During  the 
growth  of  a  GaAs  over-layer,  its  morphology  is  af¬ 
fected  by  the  interaction  of  the  inhomogeneous  strain 
in  the  neighborhood  of  the  dots,  as  well  as  in  the  wet¬ 
ting  layer.  This  interaction  affects  the  apex  of  the  dots. 
Thus,  there  is  a  thermodynamically  favored  tendency 
for  the  adatoms  at  the  apex  to  migrate  to  the  side  of 
the  dots  [15].  One  consequence  of  this  migration  is  a 
reduction  in  the  height  of  the  buried  dot  and  formation 
of  a  flat  (0  0  1 )  surface. 

As  stated  earlier,  our  nanostructure  samples  were 
doped  to  yield  about  2  electrons  per  dot.  Because 
of  tunneling  processes,  the  ground  states  of  the 
large-size  dots  are  the  ones  filled  by  most  of  the  avail¬ 
able  electrons.  The  ground  states  of  the  small-size 
dots,  on  the  other  hand,  may  be  partially  filled  or  not 
filled  at  all.  Thus,  the  intraband  absorbance  of  the 
small-size  dots  is  likely  to  be  less  than  that  of  the 
large-size  dots  at  lower  temperatures.  As  the  temper¬ 
ature  is  raised,  a  redistribution  of  the  carriers  takes 
place.  Under  equilibrium  conditions,  carrier  transfer 
from  large-  to-small-size  dots  is  the  most  favored 
process.  As  a  consequence,  the  number  of  electrons 
in  the  large-size  dots  decreases  with  increasing  tem¬ 
perature.  This,  in  turn,  reduces  the  integrated  ab¬ 
sorbance  as  observed  in  our  experiments.  Since  there 
are  now  more  electrons  in  the  small-size  dots  (as 
temperature  is  increased),  the  integrated  absorbance 
increases;  again,  in  agreement  with  our  experimental 
observations. 

The  spectral  widths  of  the  absorption  bands  shown 
in  Fig.  6  are  fairly  broad.  In  general,  the  spectral 


spreads  of  interband  and  intraband  transitions  are  a  re¬ 
sult  of  homogeneous  and  inhomogeneous  broadening 
processes.  Homogeneous  broadening,  which  is  due  to 
different  types  of  scattering  mechanisms,  usually  gives 
rise  to  a  Lorentzian  line  shape.  Inhomogeneous  broad¬ 
ening,  on  the  other  hand,  is  a  consequence  of  the  size 
distribution  in  the  dot  population.  This  type  of  broad¬ 
ening  gives  rise  to  a  Gaussian  line  shape.  The  emis¬ 
sion  and  absorption  spectra  observed  here  could  not 
be  fitted  to  a  single  type  of  line  shape,  Lorentzian  or 
Gaussian.  Only  a  mixture  of  both  line  shapes  gives 
spectra  that  agree  with  the  experimental  and  decon¬ 
volved  analytical  data.  This  fact  supports  the  no¬ 
tion  that  there  is  indeed  a  size  distribution  in  the  dot 
population. 

In  summary,  we  have  studied  the  effects  of  carrier 
redistribution  with  temperature  in  dots  of  different 
sizes  by  observing  the  manifestations  of  interband  and 
intraband  transitions  via  luminescence  and  absorption 
experiments.  These  experiments  are  important  in  the 
understanding  of  the  behavior  of  new  or  improved 
optoelectronic  devices  fabricated  from  quantum-dot 
nanostructures  intended  for  room-temperature 
operation. 
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Abstract 

We  have  optically  characterized  the  emission  of  a  layer  of  InAs  quantum  dots  embedded  at  the  center  of  symmetric 
(In,Ga)As  quantum-well  structures.  By  grading  the  profile  of  the  indium  composition  in  the  quantum-well  regions  in  a 
certain  way,  it  is  found  that  the  room-temperature  emission  of  the  dots  can  be  enhanced  by  almost  an  order  of 
magnitude  compared  to  dots  bounded  by  conventional  (In,Ga)As  quantum  wells.  The  largest  enhancement  is  obtained 
for  structures  where  the  grading  is  achieved  by  using  a  GaAs/InAs  superlattice  structure  to  form  a  quasi-linear 
distribution  of  the  indium.  The  dots  inside  the  graded  quantum-well  structures  are  engineered  to  emit  at  the  important 
telecommunication  wavelength  of  1.3  pm.  ©  2002  Published  by  Elsevier  Science  B.V. 

PACS:  73.61. Ey;  73.63.Kv;  68.35.Bs;  78.66.Fd 

Keywords;  Al.  Low  dimensional  structures;  Al.  Optical  microscopy;  A3.  Molecular  beam  epitaxy;  Bl.  Nanomaterials;  B3.  Laser 
diodes 


The  first  GaAs-based  quantum-dot  laser  was 
reported  in  1994  [1].  Since  then,  remarkable 
progress  has  been  made  in  the  development  of 
this  device.  Edge-emitting  lasers  operating  from 
1.0  to  1.3  pm  with  very  low-threshold  currents 
have  been  reported  [2,3,9];  in  addition,  vertical- 
cavity  surface-emitting  lasers  have  been  success¬ 
fully  demonstrated  [4].  One  of  the  advantages  of 
utilizing  GaAs-based  quantum  dots  for  fabricating 
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infrared  lasers  (at  1.3  pm)  is  that  the  mature 
manufacturing  and  processing  technology  for  this 
materials  system  can  be  brought  to  bear  in  the 
production  of  these  lasers  for  data  and  telecom¬ 
munication  applications. 

There  are  currently  several  approaches  utilized 
in  the  (molecular  beam  epitaxial)  synthesis  of 
(In,Ga)As/GaAs  quantum  dots  that  emit  at 
1.3  pm.  These  include  (i)  the  alternate  supply  of 
the  group-III  indium  and  the  group-V  arsenic 
source  materials  to  form  the  InAs  dots  [5],  (ii)  slow 
growth  (<0.01  ML/s)  of  the  InAs  dots  at  high 
substrate  temperatures  [6],  (iii)  burying  the  InAs 
dots  with  an  (In,Ga)As  (rather  than  a  GaAs) 
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overlayer  to  reduce  the  strain  [7],  and  (iv) 
embedding  the  InAs  dots  inside  an  (In,Ga)As 
quantum- well  structure  [8].  The  basic  physical 
effect  in  the  first  two  approaches  is  the  formation 
of  large-size  dots,  which  result  in  a  small  effective 
quantum-dot  bandgap  whose  consequence  is  long- 
wavelength  emission  when  the  dots  are  optically  or 
electrically  pumped.  In  approach  (iii)  and  (iv),  the 
(In,Ga)As  layers  play  the  role  of  reducing  the 
lattice-mismatch-induced  strain  around  the  InAs 
dots. 

Most  quantum-dot  lasers  operating  at  1.3  pm 
appear  to  have  two  major  unresolved  problems:  a 
low  density  of  dots  and  a  weak  carrier  confinement 
to  the  active  region.  Both  of  these  problems  imply 
low  optical  gain.  In  fact,  the  low  gain  is  believed  to 
prohibit  the  lasers  from  operating  at  the  ground 
state  without  high-reflectivity  coatings  [3,10].  One 
approach  to  achieving  tight  carrier  confinement  is 
to  use  growth  method  (iv),  which  embeds  the  dots 
in  an  (In,Ga)As  quantum-well  layer.  In  addition  to 
serving  as  a  strain  reducing  (and  hence  a  wave¬ 
length  tuning)  layer,  the  (InGa)As  quantum  well 
also  serves  to  confine  carriers  within  the  vicinity  of 
the  quantum-dot  layer,  thus  promoting  capture 
within  the  dots  [9], 

In  this  paper,  we  show  that  additional  deliberate 
strain,  and  hence  band  structure  engineering  of  the 
InAs/(In,Ga)As  quantum-dot  and  (In,Ga)As/ 
GaAs  quantum-well  layers  can  result  in  a  con¬ 
trollable  tuning  of  the  emission  wavelength; 
furthermore,  these  structures  promote  an  increase 
in  the  luminescence  efficiency. 

The  samples  in  our  study  were  grown  in  a 
solid-source  molecular-beam  epitaxy  system  on 
semi-insulating  (0  0 1)  GaAs  substrates.  The  basic 
structure  consists  of  InAs  quantum  dots  symme¬ 
trically  sandwiched  between  two  (In,Ga)As  layers; 
this  structure,  in  turn,  is  inserted  at  the  center  of  a 
180-nm-GaAs  layer  on  top  of  a  250-nm-GaAs 
buffer  layer  on  a  (00 1)  GaAs  substrate.  The  250- 
nm-GaAs  buffer  layer  was  grown  at  a  substrate 
temperature  of  580°C.  The  temperature  was  then 
gradually  lowered  to  490°C  [17]  during  the  growth 
of  the  90-nm-GaAs  layer  which  precedes  the 
growth  of  the  bottom  InGaAs  layer  on  top  of 
which  the  InAs  dots  are  grown.  The  rest  of  the 
structure  was  grown  at  490°C.  Three  samples, 


different  in  the  essential  details  of  their  band 
structures,  were  grown.  A  schematic  of  the  generic 
structure  is  shown  in  Fig.  1  for  all  three  samples. 
Sample  A  consists  of  2.5  ML  of  InAs  quantum 
dots  sandwiched  at  the  center  of  an  8-nm- 
In0  ]5Ga0.85As  layer.  Sample  B  consists  of  2.5  ML 
of  InAs  quantum  dots  sandwiched  at  the  center  of 
an  8-nm-InYGa]_xAs  layer;  the  bottom  half  of  the 
InAGai_xAs  layer  is  graded  from  x  =  0.03  to  0.3; 
the  top  half  is  graded  from  x  =  0.3  to  0.03.  The 
change  of  indium  composition,  x,  in  the  graded 
In^-Gai^As  quantum- well  region  is  achieved  by 
gradually  increasing  or  decreasing  the  indium  cell 
temperature  as  a  function  of  time.  The  structure  of 
sample  C  is  similar  to  that  of  sample  B  except  that 
InAs/GaAs  superlattice  layers  are  used  to  form  the 
quasi-linear  indium  distribution  along  the  growth 
direction  (see  Fig.  1).  The  InAs/GaAs  superlattice 
we  use  in  sample  C  consists  of  nine  periods  of 
InAs/GaAs,  where  the  thickness  of  the  GaAs  layer 
in  each  period  is  4  A;  the  thickness  of  the  InAs 
layer  in  each  period  is  determined  by  the  temporal 
duration  of  each  period,  which  is  increased  in  an 
ordinal  sequence  from  1  to  9  s,  or  decreased  from  9 
to  1  s.  The  growth  rate  of  the  InAs  dots  was 
0.05  ML/s,  and  the  As4  over-pressure  was  main¬ 
tained  at  4.0  x  10~6Torr  during  the  growth  of  all 
the  samples.  The  substrate  holder  was  continu¬ 
ously  rotated  to  improve  the  uniformity.  Another 
set  of  samples,  grown  under  identical  conditions, 
but  without  the  layers  that  cover  the  dots  was 
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Fig.  1.  Schematic  of  the  conduction  band  structure  of  quan¬ 
tum-dot  samples  embedded  at  the  center  of  quantum-well 
structures.  The  inset  shows  the  physical  layer  structure. 
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grown  for  atomic  force  microscope  (AFM)  surface 
studies.  We  want  to  comment  that  during  the 
epitaxy  of  samples  B  and  C,  the  reflection  high- 
energy  electron  diffraction  (RHEED)  pattern  was 
observed  to  gradually  dim  during  the  growth  of 
the  graded  In^Ga^.As  layer  adjacent  to  the  GaAs 
buffer.  For  sample  A,  however,  the  RHEED 
intensity  was  abruptly  dimmed  after  a  few  seconds. 
This  sample  consisted  of  an  Ino.15Gao.85As  layer 
grown  directly  adjacent  to  the  GaAs  buffer  layer 
without  any  intermediate  grading. 

Our  samples  were  characterized  by  photolumi¬ 
nescence  spectroscopy  at  both  room  temperature 
(295  K)  and  low  temperature  (78  K).  All  the 
samples  were  mounted  on  a  flat  copper  plate 
whose  temperature  could  be  controlled.  An  Argon 
ion  laser,  emitting  at  488  nm,  was  used  as  an 
excitation  source;  the  emitted  radiation  was 
detected  with  a  cooled  Ge  detector. 

The  room-temperature  photoluminescence  spec¬ 
tra  of  samples  A,  B,  and  C  are  shown  in  Fig.  2. 
The  emission  from  the  ground-state  transition  for 
all  three  samples  is  at  or  beyond  1.3  pm,  with 
sample  C  extending  to  1.35  pm.  The  spectral 
linewidths,  at  half-maximum  intensity,  for  the 
samples  are  Ava  =  48  meV,  Avb  =  53  meV,  and 
Avc  =  39  meV,  respectively.  These  linewidths  are 
fairly  broad,  indicating  the  effects  of  inhomoge¬ 
neous  broadening  due  to  the  variation  of  the  dot- 


Fig.  2.  Room-temperature  photoluminescence  spectra  of  the 
quantum-dot  samples  at  an  optical  excitation  level  of  ~30W/ 


sizes.  The  grading  of  the  InxGai_xAs  quantum- 
well  confining  layers  does  not  have  any  apparent 
effect  on  the  uniformity  of  the  dots-sizes.  What  the 
grading  does  appear  to  do  is  to  improve  the 
efficiency  of  emission.  From  Fig.  2,  it  is  evident 
that  there  is  an  enhancement  of  the  photolumines¬ 
cence  emission  intensity  for  samples  B  and  C.  All 
three  samples  were  excited  with  the  same  flux, 
under  the  same  testing  conditions.  The  peak 
intensity  for  sample  C,  for  example,  is  about  8.7 
times  higher  than  that  for  sample  A.  We  have  also 
performed  X-ray  diffraction  measurements  on  the 
samples  to  evaluate  the  overall  strain  in  the 
structures  using  the  [004]  reflection.  The  X-ray 
measurements  were  performed  using  a  Phillips 
MRD  system.  The  average  strain  in  the  structures, 
parallel  to  the  growth  direction,  was  determined  as 
1.12  x  1(T2  for  sample  A,  1.02  x  10-2  for  sample 
B,  and  1.13  x  10~2  for  sample  C.  The  correspond¬ 
ing  indium  compositions  were  15.7%  for  sample 
A,  14.3%  for  sample  B,  and  15.8%  for  sample  C. 
The  slightly  lower  indium  composition  for  sample 
B  is  due  to  the  non-linear  relationship  between  the 
indium  cell  temperature  and  the  flux  arrival  rate  at 
the  growth  surface  during  the  grading.  Since  in 
these  samples,  the  only  significant  difference  is  the 
structure  of  the  In*Gai_*As  layer  into  which  the 
quantum  dots  are  embedded,  it  is  reasonable  to 
speculate  that  the  enhancement  in  luminescence 
intensity  is  due  to  the  nature  of  the  grading.  There 
may  be  two  possible  reasons  that  account  for  this 
enhancement.  One  is  that,  in  a  graded  QW,  the 
energy  band  may  also  be  graded;  this  leads  to  the 
existence  of  a  field  which  helps  accelerate  the 
carriers  along  the  quantum-well  gradient,  and 
hence  into  the  quantum  dots.  Because  of  the 
three-dimensional  nature  of  quantum  dots,  the 
graded  potential  also  exists  in  the  lateral  direction; 
this  enhances  the  spatial  carrier  capture  within  the 
dots.  The  second  reason  for  the  enhancement  of 
luminescence  is  the  smoothing  effect  of  the  InAs/ 
GaAs  superlattice.  The  superlattice  causes  the 
transition  of  the  lattice  parameter  from  that  of 
GaAs  to  the  one  InAs  to  be  continuously  smooth, 
thus  reducing  the  formation  of  lattice-mismatch- 
induced  defects.  The  defects,  if  present,  would 
form  loss  channels  in  the  electron-hole  recombi¬ 
nation  process.  Of  the  two  grading  techniques 
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Fig.  3.  The  Arrhenius  plots  of  the  integrated  luminescence 
intensity  of  samples  A,  B  and  C  for  the  optical  excitation  level 
of  ~10W/cm2.  The  solid  lines  are  fitted  to  Eq.  (1)  in  the  text. 
The  inset  shows  the  ratio  of  the  integrated  luminescence 
intensity  of  sample  B  to  sample  A,  and  of  sample  C  to  sample 
A  as  a  function  of  temperature. 

used,  the  quasi-linear  grading  method  of  sample  C 
yields  the  better  results. 

We  have  further  studied  the  temperature-depen- 
dent  integrated  luminescence  intensity  of  the 
samples.  The  Arrhenius  plots  for  the  samples  are 
shown  in  Fig.  3.  The  variation  of  the  integrated 
intensity  data  with  temperature  can  be  described 
by  the  generic  empirical  relationship  [1 1] 

I(T)  = 

- — -  (1) 

1  +  Cj  exp(-£i  fkT)  +  C2  exp {-E2/kT)'  v  ’ 

where  E\  and  Ei  are  the  thermal  activation 
energies  for  loss  mechanisms  active  at  certain 
temperature  ranges,  k  is  the  Boltzmann  constant, 
T  is  the  temperature,  and  J0,  C\  and  C2  are  (fitting) 
constants.  In  a  physical  model,  these  constants 
would  take  into  account  the  recombination  rates 
and  the  geometric  dimensions  of  the  dots.  The 
fitting  parameters  and  extracted  activation  ener¬ 
gies  are  given  in  Table  1.  In  Fig.  3,  the  solid  lines 
show  the  fit  to  the  data  using  Eq.  (1)  given  above; 
the  filled  circles,  open  circles,  and  squares  repre¬ 
sent  experimental  data.  The  extracted  activation 
energies  are  E\  =  289  meV,  and  £2  =  41  meV  for 
sample  A;  E\  =  263  meV,  and  £2  =  33  meV  for 
sample  B;  and  E\  —  255  meV  and  £2  =  25meV 


Table  1 


Sample 

Fitting  parameters 

Activation  energies 
(meV) 

h 

x  ~ 

o 

C2 

Ex 

e2 

A 

0.252 

14.8 

32.7 

289 

41 

B 

0.334 

1.34 

18.6 

263 

33 

C 

0.532 

0.706 

9.31 

255 

25 

for  sample  C.  Note  that  the  luminescence  intensity 
has  a  quenching  threshold  temperature;  for  both 
sample  B  and  C  this  temperature  is  — 190  K,  and  it 
is  ~170K  for  sample  A.  The  loss  channel  that 
dominates  at  these  high  temperatures  is  probably 
associated  with  the  activation  energy  E\.  A 
calculation  (which  we  will  discuss  later)  shows 
that  the  energy  difference  between  the  ground  state 
in  the  conduction  band  of  a  dot  and  the  band-edge 
of  the  surrounding  (In,Ga)As  quantum  well  at 
300  K  is  about  317  meV  for  structure  A  and 
267  meV  for  sample  C.  These  energy  differences 
are  quite  close  to  the  measured  activation  energies 
E\  for  samples  A  and  C.  At  the  quenching 
threshold  temperatures,  the  thermal  escape  of 
carriers  from  the  dots  begins  to  play  an  important 
role  in  carrier  recombination  dynamics.  It  is  then 
conceivable  that  the  carriers  have  sufficient  energy 
to  overcome  these  barriers.  An  interesting  aspect 
of  the  comparison  between  the  recombination 
dynamics  of  our  samples  is  that  the  photolumines¬ 
cence  enhancements  of  samples  C  and  B  are  more 
dramatic  at  the  higher  temperatures  (see  the  inset 
in  Fig.  3).  Our  results  here  are  similar  to  those  of 
Ru  et  al.  who  used  hydrogen  passivation  (in  order 
to  reduce  the  defect  concentrations  in  their 
samples)  to  enhance  the  photoluminescence  emis¬ 
sion  [12],  In  our  case,  grading  of  the  QW  (and  the 
use  of  the  superlattice  in  sample  C)  may  reduce  the 
tendency  for  lattice-mismatch-induced  formation 
of  defects  at  the  QW/QD  interface  region,  thus 
reducing  non-radiative  defect-related  losses  in  the 
structures. 

We  have  examined  the  possible  effects  of  the 
grading  on  the  dot  density  and  uniformity.  Three 
structures  similar  to  samples  A,  B  and  C  were 
grown  for  AFM  observations.  Fig.  4  shows  the 
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(a)  Sample  A,  1 .25  by  1 .25  (b)  Sample  B,  1.25  by  1.25  jim 


(c)  Sample  C,  1.25  by  1.25  pm 


Fig.  4.  AFM  photo-micrographs  of  the  morphology  of  samples  A,  B,  and  C. 


AFM  micrographs  of  the  three  samples.  The  areal 
dot  densities  of  these  samples  are  2.5  x  10locm~2 
for  sample  A,  1.8  x  10locm-2  for  sample  B,  and 
2.8  x  10locm-2  for  sample  C.  The  true  size  of 
the  dots  cannot  be  determined  from  the  AFM 
measurements;  reasonable  estimates,  however,  can 
be  obtained  from  cross-sectional  TEM  studies. 

To  show  the  basic  quantum-confined  features  of 
the  band  structure  of  the  samples,  we  have  studied 
the  intensity-dependent  photoluminescence  emis¬ 
sion  at  both  room  and  LN2  temperatures.  In 
Fig.  5,  we  show  the  LN2  (78  K)  photoluminescence 
emission  for  the  three  samples  at  high  pump  levels. 
There  are  clearly  four  transition  energy  peaks  for 
each  sample.  We  have  carried  out  some  analytical 
calculations  to  determine  the  locations  of  the 
interband  transition  energies  for  samples  A  and 
C.  Because  of  the  non-linear  composition  of  the 
indium  in  the  QW  region  of  sample  B,  we  have,  for 
simplicity,  omitted  it  in  our  calculations.  These 


Fig.  5.  The  low-temperature  (78  K)  photoluminescence  spectra 
of  samples  A,  B,  and  C  under  high  optical  excitation. 


calculations  require  information  on  the  strain 
distribution  inside  and  around  the  dots.  The  strain 
generally  depends  on  the  shape  and  size  of  the 
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Table  2 


Structure  Data  source  Transition  energy  (eV) 


Low  temperature  (77  K) 

RT  (300  K) 

Eo 

Ei 

Ei 

Ei 

Eo 

A 

Experimental 

1.019 

■m 

■m 

0.934 

Calculated 

1.0162 

iUfa 

0.9405 

C 

Experimental 

0.987 

1.036 

1.083 

1.139 

0.919 

Calculated 

0.9830 

1.0379 

1.0809 

1.1315 

0.9138 

dots,  which  are  difficult  to  determine  accurately  in 
most  experimental  situations.  For  our  dots,  we 
estimate  their  lateral  extent  (base)  to  be  around 
15nm;  the  base-to-height  ratio  for  sample  A  is 
about  7;  and  it  is  5  for  sample  C.  These  estimates 
are  based  on  our  earlier  work  on  buried  InAs / 
GaAs  and  InGaAs/GaAs  quantum-dot  work  [13], 
and  on  TEM  and  SEM  observations.  The  estimate 
on  height  is  further  corroborated  with  in-situ 
RHEED  observations  where  the  thickness  of  an 
over-layer  needed  to  completely  bury  an  array  of 
dots  (so  that  the  layer-by-layer  growth  mode  is  re¬ 
established)  can  give  a  rough  measure  of  the  dot 
height.  For  the  calculations,  we  have  further 
assumed  that  the  dots  are  truncated  pyramids. 

The  electronic  spectra  here  are  calculated  in  the 
envelope  function  approximation  using  an  eight- 
band,  strain-dependent  Hamiltonian  based  on  the 
k « p  method  [14].  The  details  of  the  method  we  use 
in  the  calculations  have  been  published  elsewhere 
[15,16].  The  calculated  transition  energies  for  our 
samples  A  and  C  are  shown  in  Table  2.  We  have 
also  shown  the  experimentally  determined  values 
from  the  photoluminescence  data  obtained  from 
multi-peak  Gaussian  fitting.  The  experimental  and 
the  calculated  transition  energies  are  in  reasonably 
good  agreement. 

In  summary,  we  have  shown  that  the  photo¬ 
luminescence  efficiency  of  (In,Ga)As/GaAs  quan¬ 
tum  dots  can  be  enhanced  by  embedding  them  in  a 
graded  (In,Ga)As  quantum-well  structures.  We 
believe  this  approach  is  simpler  and  more  practical 
than  the  use  of  in  situ  hydrogen  passivation  for 
defect  annealing.  Furthermore,  our  proposed 
approach  has  the  virtue  of  promoting  carrier 


capture  in  potential  laser  structures  in  addition 
to  allowing  emission  at  the  desirable  telecommu¬ 
nications  wavelength  of  1.3  pm.  The  band  grading 
we  describe  here  is  reminiscent  of  the  graded  index, 
separate  confinement  heterostructures 
(GRINSCH)  of  early  generation  active  regions  of 
quantum-well  lasers. 

This  work  is  supported  under  Grant  Number 
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This  paper  presents  a  study  of  electronic  carrier  redistribution  in  doped  (In,  Ga)As/GaAs  quantum  dots  due  to  size  variation. 
The  dots  in  our  samples  fall  roughly  into  “large”  and  “small”  size  clusters.  At  low  temperatures,  the  carriers  predominantly 
occupy  the  large- size  dots.  As  the  sample  temperature  is  raised,  thermal  energy  redistributes  the  carriers  so  that  some  begin  to 
occupy  energy  states  in  the  small-size  dots.  At  room-temperature,  the  photoluminescence  emission  spectrum  exhibits  two 
distinct  peaks  due  to  transitions  from  the  large-  and  small-size  dots.  This  carrier’ redistribution  can  also  be  followed  by 
Observing  the  integrated  absorbance  as  a  function  of  temperature.  The  transition  energy  levels  extracted  from 
photoluminescence  data  compare  favorably  with  theoretically  calculated  values  for  dots  with  truncated  pyramidal 
shapes.  [DOI:  10.1 143/JJAP.4 1.482] 

KEYWORDS:  quantum  dots,  intraband  absorption,  photoluminescence,  (In,  Ga)As/GaAs,  nanostructures,  strain,  carrier 
redistribution,  quantum  dot-size  distribution 


1.  Introduction 

The  study  of  semiconductor  quantum-dot  nanostructures 
has  attracted  considerable  interest  in  the  last  few  years.  This 
is  because  of  the  potential  application  of  the  dots  in 
electronic  and  optoelectronic  devices.  These  nanostructures 
confine  carriers  in  all  three  dimensions.  As  a  consequence  of 
this,  several  important  characteristics  emerge:  these  include 
atom-like  density  of  states,  large  exciton  binding  energies, 
and  enhanced  oscillator  strengths.  Some  of  these  character¬ 
istics  have  enabled  the  demonstration  of  novel  light-emitting 
sources,^  single-electron  memory  elements, 2)  and  infrared 
photodetectors.3)  Of  the  devices  demonstrated  so  far, 
perhaps  the  most  advanced  is  the  interband  laser; 1,4)  it  has 
operated  at  both  low-temperatures  and  room-temperature. 
The  interband  quantum-dot  laser  has  also  been  shown  to 
exhibit  a  high  characteristic  temperature,  To — which  is  a 
measure  of  the  sensitivity  of  the  threshold  current  to 
temperature  variations. 

The  intraband  absorption  properties  of  the  dots  have  been 
used  to  design  and  fabricate  infrared  photodetectors.  These 
devices,  however,  have  only  operated  at  cryogenic  tempera¬ 
tures  so  far.5,6)  Further  development  for  high-temperature 
operation  will  require  an  understanding  of  the  effects  of 
temperature  on  the  intraband  absorption  characteristics.  In 
fact,  a  better  understanding  of  the  effects  of  temperature  on 
both  the  interband  and  intraband  emission  and  absorption 
characteristics  of  the  dots  would  be  helpful  in  the  design  of 
future  room-temperature  devices.  Several  workers  have 
reported  the  effects  of  temperature  on  interband  emission 
from  InAs/GaAs,  (In,  Ga)As/GaAs,  and  (Al,  In)As/(Al,  - 
Ga)As  quantum-dot  structures  .7-9)  It  has  been  found  that  in 
some  cases,  the  emission  intensity  from  the  nanostructures  is 
quenched  as  the  temperature  is  raised.  When  dots  of 
different  sizes  exist  in  an  ensemble,  the  variation  in 
temperature  has  been  linked  to  a  thermal  redistribution  of 
carriers  within  the  ensemble.  These  observations  were  made 
on  undoped  quantum  dots;  no  studies  on  doped  nanostruc- 
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tures  have  yet  been  reported.  The  emission  and  absorption 
behavior  of  doped  quantum  dots,  as  the  temperature  is 
changed,  is  important  because  it  has  a  direct  bearing  on  the 
operating  characteristics  of  many  optoelectronic  devices. 

While  most  of  the  work  in  the  literature  on  quantum  dots 
has  focused  on  dot  emission  characteristics,  there  is  now  a 
growing  interest  in  absorption:  the  inverse  of  emission. 
Absorption  provides  a  unique  signature  for  the  energy 
structure  in  the  quantum  dots.  Normal-incidence  intraband 
absorption,  for  example,  can  provide  information  on 
intersublevel  energy  separations.  Several  groups  have 
reported  low-temperature  photoconduction  spectroscopy 
measurements  in  n-type,  doped  (In,  Ga)As/GaAs3)  and 
InAs/GaAs 10)  dots.  The  work  on  InAs  dots  has  been 
performed  at  both  low-  and  room-temperatures. n,12)  Since 
most  dot  systems  studied  to  date  contain  distinct  size 
distributions,  it  would  be  interesting  to  know  the  effect  of 
thermal  carrier  redistribution  on  the  absorption  strengths  of 
the  dot  systems. 

One  of  the  most  challenging  aspects  of  quantum-dot  work 
is  the  development  of  quantitative  models  that  accurately 
describe  the  dots  synthesized  from  m-V  compound 
semiconductors.  A  number  of  elegant  theoretical  models 
exist; 13,14)  the  shapes  of  the  dots  and  the  parameters  used  in 
the  theoretical  calculations,  however,  often  yield  results  that 
are  difficult  to  match  to  experimental  data.  As  a  conse¬ 
quence,  the  models  have  not  played  a  significant  role  in 
guiding  experimentation.  Several  groups  are  therefore 
engaged  in  the  development  of  models  that  best  suit  their 
experimental  situation.  In  this  work,  we  use  a  transmission 
electron  microscope  (TEM)  to  experimentally  determine  the 
shape  of  the  dots.  Our  observations  lead  us  to  believe  that 
the  dots  grown  in  our  laboratory  are  predominantly  truncated 
pyramids.  We  use  the  dimensions  and  shape  extracted  from 
the  TEM  images  in  calculations  of  the  energy  levels  in  the 
dots.  In  these  calculations,  one  must  take  into  account  the 
effects  of  strain — which,  in  the  first  place,  is  responsible  for 
the  formation  of  the  dots.  Two  of  the  most  common  methods 
used  in  the  quantitative  analysis  of  strain  involve  spatial 
discretization  of  the  region  within  which  the  strain  is  known 
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to  exist:  the  first  method  is  the  finite  element  technique14^ 
and  the  second  is  the  finite  difference  method. 15)  Both 
techniques  require  considerable  computational  resources.  In 
this  work  we  will  use  a  simple  analytical  technique  for 
calculating  the  strain  distribution  in  and  around  the  dot.  The 
results  of  the  strain  calculations  are  a  precursor  to  the  energy 
structure  analysis. 

Our  aim  in  this  paper  is  to:  (a)  show  the  effects  of 
temperature  on  the  luminescence  characteristics  of  Si-doped 
(In,  Ga)As  quantum  dots  and  the  associated  Boltzmann 
redistribution  of  carriers  in  dot  ensembles  with  size 
variations  in  them;  (b)  show  that  the  doped  dots  exhibit 
strong,  normal -incidence,  intraband  absorption;  and  (c) 
compare  the  experimentally  determined  energy  levels  in 
our  quantum  dots  with  values  derived  from  calculations. 

2.  Experimental 

The  quantum- dot  superlattices  used  in  this  study  were 
grown  by  a  solid-source  molecular  beam  epitaxy  system.  A 
typical  structure  consists  of  five  periods  of  (In,  Ga)As/GaAs 
quantum  dots  grown  on  a  (001)  GaAs  3ub3trntc.  The  nominal 
indium  composition,  based  on  growth  parameters  and  in-situ 
RHEED  observations,  was  about  30%.  X-ray  analysis  of  the 
samples,  however,  indicated  that  the  average  indium 
composition  in  the  dots  was  about  23%.  The  dots  were  8- 
doped  with  silicon  to  yield  about  two  electrons  per  dot  The 
thickness  of  the  GaAs  barrier  layer  in  the  superlattice  was 
about  50  nm.  An  uncapped  sample  was  also  grown  using 
identical  growth  parameters  for  surface  analysis  of  the  dots 
by  an  atomic  force  microscope.  The  surface  studies  indicate 
that  the  density  of  the  dots  is  about  6.3-6. 8  x  1010  cm"2,  and 
from  high-resolution  transmission  electron  microscopy,  we 
were  able  to  determine  that  the  buried  dots  are  about  3  nm 
high;  the  lateral  extent  of  their  bases  ranged  from  about  16  to 
20  nm. 

2.1  Photoluminescence 

Photoluminescence  measurements  were  carried  out  using 
an  argon-ion  laser  ( X  —  488  nm)  as  the  excitation  source. 
The  temperature  of  the  samples  was  varied  from  77  to  300  K. 
The  emitted  radiation  was  detected  with  a  cooled  Ge 
detector.  Figure  1  shows  the  luminescence  spectra  obtained 
as  the  temperature  is  varied  from  liquid-nitrogen  (77  K)  to 
room-temperature  (300  K).  Note  the  quenching  of  the 
peak  intensity  as  the  temperature  is  raised;  also  observe 
that  as  the  temperature  approaches  300  K,  the  spectrum 
broadens  and  there  is  a  suggestion  of  a  second  peak  hidden 
by  the  broad  emission.  This  is  clearly  shown  in  the  room- 
temperature  spectrum  of  Fig.  2;  at  low-temperature  (77  K), 
however,  the  emission  is  characterized  by  a  single  peak 
located  at  1.19  eV,  with  a  spectral  spread  of  about  30.8  meV. 
Furthermore,  no  luminescence  from  the  states  in  the  wetting 
layer  was  observed  due  to  photo-generated  carriers  in  the 
wetting  layer  transfer  into  the  adjacent  dots.  A  deconvolu¬ 
tion  of  the  room-temperature  spectrum  indicates  that  it  is, 
indeed,  a  composite  of  two  spectra  with  peaks  located  at 
1.12  and  1.16  eV.  These  spectra  have  been  labeled  (a)  and 
(b)  in  Fig.  2.  The  energy  difference  between  the  peaks  is 
about  40meV;  this  is  much  smaller  than  the  nominal 
separation  (80-100  me V)  between  the  ground  state  and  the 
first  excited  state  of  the  conduction  band  of  our  quantum 


dots.  From  intraband  absorption  measurements  (to  be 
discussed  later),  the  nominal  separation  should  be  about 
100  meV.  The  appearance  of  the  second  peak — on  the  high 
energy  side  of  the  room-temperature  spectrum — therefore 
cannot  be  attributed  to  an  excited  state  transition.  One  can 
explain  the  appearance  of  the  second  peak  by  considering  a 
bimodal  size  distribution  for  the  dots.16)  On  this  premise,  we 
assume  that  the  dots  either  belong  to  a  “large-size”  group  or 
a  “small-size”  group.  The  origin  of  the  two  peaks  is 
therefore  a  consequence  of  thermal  energy  redistribution  of 
carriers  among  the  dots  of  the  two  groups.  The  qualitative 
descriptors  of  large-  and  small-size  represent  average  dot 
dimensions.  At  low  temperatures  and  moderate  excitations, 
the  source  generates  electron-hole  pairs  that  predominantly 
tunnel  to  occupy  the  ground  states  of  the  “large-size”  dots 
since  these  represent  energy  minima  for  the  ensemble.  This 
would  explain  why  at  77  K,  only  a  single  emission  peak  is 


Fig.  1.  Photoluminescence  spectra  of  (In,  Ga)As/GaAs  quantum  dots  at 
different  sample  temperatures. 


Fig.  2.  Luminescence  spectra  of  (In,  Ga)As/GaAs  quantum  dots  at  77  and 
300  K. 
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observed  at  1.19  eV.  As  the  temperature  is  raised,  however, 
thermal  energy  redistributes  the  carriers  among  the  dots  of 
different  sizes,  causing  some  to  occupy  the  slightly  higher 
ground  states  of  the  “small-size”  dots.  This  becomes  evident 
at  room-temperature  where  the  emission  now  is  from  both 
the  “large-size”  and  “small-size”  dots,  leading  to  peaks  at 
1.12  and  1.16  eV.  We  interpret  the  higher  energy  emission  to 
be  from  the  “small-size”  dots  where  the  conduction 
(valence)  band  ground  state  energies,  on  the  average,  are 
positioned  higher  (lower)  than  those  in  the  “large-size”  dots. 

The  thermal  redistribution  of  carriers  in  the  dot  ensembles 
can  also  be  followed  in  another  way.  One  can  relate  a 
relative  shift  of  emission  energies,  A£qD(T),  from  a  fixed 
reference  point  to  a  variation  in  temperature.  We  consider  a 
single  low-temperature  emission  peak  fcv(107K)  as  the 
reference  point.  The  relative  shift  is  therefore  defined  as: 

AEqd(T)  =  hv(T)  -  hv{  107  K),  (1) 

where  hv(T)  is  the  temperature-dependent  peak.  The  relative 
shifts  of  the  two  spectra  (a)  and  (b)  of  Fig.  2,  as  functions  of 
temperature,  are  shown  in  Fig.  3.  For  comparison,  we  also 
show,  as  a  solid  line,  the  relative  shift  of  the  band  gap  energy 
of  the  Ino.23Gao.77  As  material  out  of  which  the  quantum  dots 
are  synthesized.  This  variation  is  calculated  from 

AEg(T)  =  Eg(T)-Eg(101K ),  (2) 

with  the  temperature  variation  of  the  band  gap,  Eg(T ), 
calculated  according  to  Varshini’s  law.17)  In  general,  the 
peak  luminescence  energy  of  the  quantum-dot  ensemble 
should  follow  the  thermal  shrinkage  of  the  band  gap  as  the 
temperature  is  raised,  as  long  as  there  is  no  transfer  of 
carriers  from  one  state  to  another.  However,  from  Fig.  3,  one 
observes  that  the  relative  rate  of  thermal  shrinkage  for 
spectra  (a)  and  (b)  is  lower  than  that  of  the  band  gap, 
indicating  a  carrier  transfer  from  low  to  high  energy  states.  If 
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Fig.  3.  Variation  of  the  relative  shift  A£qd  (from  the  fixed  reference  point 
of  the  peak  emission  at  107  K)  as  a  function  of  temperature  for  emission 
band  (a)  and  (b)  of  Fig.  2.  The  solid  line  is  the  relative  band  gap  variation 
of  (In,  Ga)As,  A£g,  with  temperature.  The  Stokes  shift  (SS),  defined  as 
the  variation  of  the  difference  between  AEqd  and  AES  with  temperature, 
is  also  shown  on  the  right  hand  side  of  the  figure. 


D.  Pal  et  al 

the  carrier  transfer  had  been  from  high  to  low  energy  states, 
the  relative  thermal  shrinkage  rate  would  have  been  larger 
than  that  of  the  band  gap.18) 

The  difference  between  AEqq  and  A Eg  can  generally  be 
considered  to  be  the  Stokes  shift  (SS)  between  the 
luminescence  and  absorption  lines.  This  shift  is  shown  in 
Fig.  3.  Note  that  the  Stokes  shift  for  spectrum  (b)  of  Fig.  2  is 
higher  than  that  for  spectrum  (a),  suggesting  that  the  size 
distribution  of  dots  corresponding  to  emission  band  (b)  is 
larger  than  that  for  band  (a).  The  inhomogeneous  broadening 
of  the  luminescence  spectrum  depends  on  the  size  distribu¬ 
tion;  a  large-size  distribution  gives  rise  to  a  broad  spectral 
spread.  This  is  evident  in  the  luminescence  spectra,  where 
the  spectral  spread  for  band  (b)  is  61.7  meV — which  is  about 
twice  the  value  (32.8  meV)  for  band  (a). 

To  gain  further  insight  into  the  thermal  energy  redistribu¬ 
tion  of  carriers  in  the  dot  ensemble,  we  have  plotted  the 
ratio,  /?,  of  the  integrated  luminescence  intensity  for  band  (a) 
and  (b)  of  Fig.  2  as  a  function  of  temperature.  This  ratio  is 
shown  in  Fig.  4.  The  solid  line  through  the  points  is  an 
Arrhenius  plot  of  the  equation 

/  A£\ 

s_Aexp(-_j.  (3) 

The  meaning  of  the  symbols  is  discussed  in  the  following. 
The  activation  energy  A E  ~  38  meV.  The  data  displayed  in 
Fig.  4  can  be  explained  if  one  assumes  a  quasi-equilibrium 
distribution  of  carriers  between  the  “large”  and  “small” 
dots  at  high  temperatures.  The  emission  intensity  from  the 
large  and  small  dots  can  be  written,  respectively,  as: 

/jarge  =  ^large^largePlarge  >  (4) 

and 

Ismail  ”  ^small wsmaltP small »  (5) 

where  n  and  p  are  the  electron  and  hole  populations  of  the 
states  in  the  dots,  and  B  is  a  constant.  If  the  carrier 
populations  obey  Boltzmann  statistics  at  any  temperature  T , 
then  the  ratio  R  can  be  written  as: 
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Fig.  4.  The  Arrhenius  plot  of  the  ratio  of  the  integrated  photolumines¬ 
cence  intensity  of  spectra  (a)  and  (b);  the  solid  line  is  a  fit  to  the  equation 
R  =  A  zxp(—AEfKvT). 
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R  = 


RnRv  exp 


(- 


AEn  +  A£p\ 
k^T  ) 


where  Rn  and  Rv  are,  respectively,  the  ratios  of  the  electron 
and  hole  effective  densities  of  states  in  the  large-  and  small- 
size  dots.  The  energy  differences  between  electron  and  hole 
levels  in  the  large-  and  small-size  dots  are,  respectively,  A En 
and  AEp.  The  energy  difference  between  the  emission  bands 
from  the  large-  and  small-size  dots  is  given  by  A E  =  A En+ 
AEP.  If  the  temperature  dependence  of  B\^ ge,  £smaii,  Rn,  and 
Rp  is  less  than  the  exponential  term,  then  the  ratio  R  can  be 
approximated  as  exp(-A E/k^T),  which  is  consistent  with 
the  results  displayed  in  Fig.  4.  Note  that  the  energy 
difference  of  40  meV  measured  between  the  peaks  of  spectra 
(a)  and  (b)  of  Fig.  2  is  quite  close  to  the  energy  (38  meV) 
extracted  as  the  activation  from  Fig.  4.  These  results  suggest 
that  the  relative  populations  of  carriers  in  the  states  of  the 
dots  of  various  sizes  obey  Boltzmann  statistics.  Similar 
results  have  been  reported  for  undoped  quantum  dot 
samples. 19) 


2.2  Intraband  absorption 

The  energy  spacing  between  sublevels  within  the 
conduction  band  of  the  quantum-dot  ensemble  can  be 
probed  by  Fourier  transform  infrared  spectroscopy.  We  have 
carried  out  such  measurements  for  our  samples  and  the  low- 
and  room-temperature  absorption  spectra  are  shown  in 
Fig.  5.  The  spectra  have  been  normalized  by  using  the 
absorption  of  a  (semi-insulating)  GaAs  substrate  without 
quantum-dots  polished  in  a  manner  similar  to  that  of  the 
quantum-dot  samples.  The  low-temperature  (77  K),  normal- 


Wavelength  (jim) 


Energy  (meV) 


Fig.  5.  Normal-incidence  intraband  absorption  spectra  of  (In,Ga)As 
quantum  dots  at  77  and  300  K.  The  spectra  (a),  (b)  and  (c),  (d)  were 
obtained  by  deconvolution  of  the  absorption  spectra  at  77  K  and  300  K, 
respectively.  The  inset  of  this  figure  shows  simplified  band  diagrams  for 
large-  and  small-size  dots.  The  energy  levels  derived  from  interband 
emission  and  intraband  absorption  measurements  are  also  shown. 


incidence  absorbance  for  a  five-period  superlattice  of 
Ino.23Gao.77As/GaAs  quantum  dots  was  about  7%.  The 
absorbance  per  quantum-dot  superlattice  period  is  therefore 
~1.4%.  Weber  et  al 12)  have  reported  a  low-temperature 
(77  K)  absorbance  per  InAs/InAlAs  quantum-dot  period  of 
0.78%.  In  another  report  Phillips  et  u/.,n)  measured  a  room- 
temperature  (300  K)  absorbance  of  about  0.48%  per  InAs / 
GaAs  quantum-dot  period.  This  absorbance  is  less  than  the 
value  of  1.02%  we  obtain  for  our  (In,  Ga)As/GaAs  quantum- 
dot  structures  at  300  K.  We  must  emphasize  that  it  is  difficult 
to  compare  absorbances  of  samples  with  different  constitu¬ 
ents  in  the  superlattice  periods. 

The  spectra  displayed  in  Fig.  5  suggest  that  the  room-  and 
low-temperature  absorption  spectra  can  be  resolved  into  two 
resonances.  At  liquid  nitrogen  temperature  (77  K),  the 
absorption  has  resonances  at  105  meV  [curve  (a)]  and 
137meV  [curve  (b)].  At  300  K,  the  absorption  curve  can 
be  resolved  to  reveal  resonances  at  107  and  135  meV.  The 
integrated  absorbance  for  the  peak  which  moves  from  105  to 
107meV  decreases  while  that  for  the  peak  that  goes  from 
135  to  137  meV  increases  as  the  temperature  is  increased 
from  77  to  300  K.  This  is  likely  to  be  due  to  the  escape  of 
carriers  from  the  ground  states  of  the  large-size  dots  to  the 
higher  lying  ground  states  of  the  small-size  dots.  The 
integrated  absorbance  of  the  large-size  dots  (at  105  and 
107meV),  therefore  decreases  as  carriers  escape  from  these 
dots  into  the  small-size  dots.  On  the  other  hand,  the 
integrated  absorbance  of  the  small-size  dots  (at  135  and 
137meV)  increases.  On  the  basis  of  our  photoluminescence 
emission  and  intraband  absorption  data,  a  simplified  energy 
structure  for  the  large  and  small  dots  can  be  derived;  such  an 
energy  scheme  is  shown  in  the  inset  of  Fig.  5.  Here,  we  have 
assumed  that  the  energy  difference  between  the  ground 
states  in  the  valence  band  of  the  large  and  small  size  dots  is 
about  2meV.  From  an  analysis  of  the  photoluminescence 
emission  intensity  as  a  function  of  temperature,  it  was  shown 
that  in  the  conduction  band,  the  ground  state  of  the  small- 
size  dots  lies  about  38meV  above  the  ground  state  of  the 
large- size  dots. 

The  spectral  spreads  of  the  absorption  spectra  in  Figs. 
5(a),  5(b),  5(c),  and  5(d)  are,  respectively,  42.5,  31.3,  54.8, 
and  41.8  meV.  These  linewidths,  like  those  for  the  interband 
transitions,  are  a  result  of  homogeneous  and  inhomogeneous 
broadening  processes  in  the  dot  ensembles.  Homogeneous 
broadening  generally  gives  rise  to  a  Lorentzian  line-shape. 
Inhomogeneous  broadening,  which  is  mostly  due  to  size 
distribution  effects  in  dot  population,  gives  rise  to  a 
Gaussian  line-shape.  The  emission  and  absorption  spectra 
obtained  for  our  quantum-dot  structures  could  not  be  fitted 
very  well  to  either  a  pure  analytic  Lorentzian  line-shape  or  a 
pure  Gaussian  line-shape.  The  experimental  spectra  could 
only  be  fitted  to  a  mixture  of  both  a  Lorentzian  and  a 
Gaussian.  This  fact  is  further  evidence  of  the  presence  of 
size  distribution  in  the  dot  ensemble. 

3.  The  Electronic  Energy  Structure  of  the  Dots 

The  energy  structure  of  the  dots  depends  on  the  strain 
distribution  around  and  inside  the  dots.  Among  the  factors 
that  affect  this  are  the  shape  and  size  of  the  dots.  The  dots 
synthesized  in  our  laboratory  have  truncated  pyramidal 
shapes.  They  typically  lie  on  top  of  a  wetting  layer  whose 
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thickness  in  the  structures  studied  here  is  about  5  nm.  The 
average  lateral  extent  of  the  dots  is  about  16-20nm;  the 
height  is  about  3  nm. 

The  first  step  in  any  theoretical  analysis  of  the  electronic 
structure  of  the  dots  is  the  determination  of  the  strain 
distribution.  The  details  of  the  analytical  approach  for 
determining  the  strain  field  of  pyramidal  quantum  dots  with 
truncated  tops  have  been  presented  elsewhere20)  and  only  the 
salient  points  will  be  given  here.  We  use  an  analytical 
approach  based  on  the  continuum  elasticity  method 
pioneered  by  Eshelby21)  for  calculating  the  strain  distribu¬ 
tion  due  to  a  uniformly  lattice-mismatched  inclusion.  First,  a 
set  of  vectors  A  as  identified  by  Downes  et  al22)  is  defined; 
this  set  of  vectors  is  defined  such  that  V-A  yields  the 
Green’s  function  for  the  stress  components,  or*ph,  corre¬ 
sponding  to  a  spherical  point  inclusion.  Thus 

cr?Ph  =  V  •  Ay.  (7) 


The  potential  components  can  then  be  written,  according  to 
Grundmann  et  al.l4)  as: 


Aijix^z)  =- 


1  Esq  X&i 

4n  1  -  v  (x2  +  y2  +  z2)3/2 


Aij(x,y,z)  =- 


1  1  Esq  xtfj  +  xjei 

2  Ait  1  —  v  (x2  +  y2  +  z2)3/2 


(8) 


where  (i,  f)  are  replaced  by  (x,  y,  z),  and  e,-  is  the  unit  vector 
in  the  i-th  direction.  Because  of  the  linear  superposition  of 
stresses, 


Oy(x,y,z)  = j I j cr*f\x  -  x0,y -y0,z-  Zo)dV(x0,y0,zo) 


-Sit 


V •  A(x  -x0,y-  y0, z  -  Zo)dV(x0, yo, Zo). 


(9) 


Applying  the  divergence  theorem,  with  (x0,  y0,  zo)  a  point 
within  the  volume  of  the  dot,  the  stress  field  due  to  an 
arbitrarily  shaped,  uniformly  lattice-mismatched  quantum 
dot  contained  within  an  isotropic,  elastic  medium  is  obtained 
by  calculating  the  following  integral,  where  (x0,  yo,  zo)  is 
now  a  point  on  the  surface  of  the  quantum  dot.  Thus 

Oy  (*, y,  z)  =  Jj  Aij(x  -x0,y-  yo, t  -  zo)dS(x0 , y0,  Zo) 
s 


-b 

~b 


(H-z) 
2 H 


<  x  < 


b 


(H-z) 
2  H 


- <  v  <  b - , 

2H  2 H 

0  <  z  <  Ht 


(11) 


where  H  is  the  height  of  the  pyramid  in  the  absence  of 
truncation,  b  is  the  base  of  the  pyramid,  t  represents  the 
degree  of  truncation  (approximately  0.64  in  the  present 
case),  and  h  =  Ht  is  the  height  of  the  truncated  pyramid.  The 
origin  of  coordinates  is  at  the  center  of  the  base  and  the  z- 
axis  is  the  growth  direction.  The  strains  follow  from 
substitution  of  stresses  into  Hooke’s  law: 

Bij  =  “  [(1  +  v)0ij  ~  SyVCTkkl  (12) 

Cj 

We  can  define  the  hydrostatic  and  biaxial  strains  as 

eh  =  Bxx  +  £yy  -j-  (13) 

and 

£b  —  Ejz  2  (£xx  £yy)-  (14) 

The  calculated  strain  tensor  for  quantum  dots  with  a  base 
width  of  16.8  and  19.7  nm,  respectively  are  shown  in  Figs.  6 
and  7,  in  the  x~z  plane.  The  effect  of  the  strain  profile  on  the 
electronic  structure  of  the  quantum  dots  can  be  qualitatively 
understood  by  examining  the  strain-modified  band  offsets. 
The  strain-induced  shift  in  the  conduction  band  depends  on 
the  hydrostatic  component  of  the  strain  through 

SEC  =  ClciSxx  +  £yy  +  %),  (15) 

where  ac  is  the  deformation  potential  for  the  conduction 
band.  The  offset  of  the  unstrained  conduction  band  as  a 
function  of  indium  mole  fraction,  x,  is  given  by,23) 

£®(x)=£jav(x)  +  ^  +  £gW  (16) 

where  A0(x)  is  the  spin-orbit  splitting,  and  Eg(x)  is  the 
unstrained  energy  band  gap.  The  conduction  and  valence 
band-edges  are  determined,  respectively,  from  the  following 
relations: 

Ec(x)  =  E?c(x)  +  8Ec(x),  (17) 

Ef(x)  =  £*av(x)  +  +  SEV, h(x)  -  l- S£v,b(x),  (18) 


+  SijYZ ~ JJf s(x~x0’y -yo,z-zo)  (io) 

V 

X  ^F(x0,yo,Zo), 

where  e0  is  the  lattice-mismatch,  E  is  Young’s  modulus,  v  is 
Poisson’s  ratio,  V  is  the  volume  of  the  quantum  dot,  and  Sy  is 
the  Kronecker  delta.  The  last  part  of  eq.  (10)  comes  from 
evaluating  the  limit  of  the  surface  integral,  as  the  field  point 
(x,  y,  z)  approaches  the  boundary  at  position  (x0,  yo,  Zo)-  The 
misfit  strain  (lattice-mismatch)  is  taken  as  negative  for  a 
material  under  compression,  such  as  InAs  on  top  of  GaAs  or 
ihjGai-jAs  on  GaAs.  The  volume  of  a  square-based, 
truncated  pyramid  is  defined  by 


and 


Ef(x)  =  £?,av(x)  -  ^  +  5£v,h(x)  -  A  SEv,b(x) 


+ 


{fhM 


-(19) 


+  A0(x)5Ev,b(x)-F--[^v,b(x)]2. 
4 


In  the  equations  above, 

^V,hW  =  civiexx  +  £yy  +  ea), 

and 


(20) 


^V,bW  —  2 bie^  -£xx~  Syy).  (21) 

£^  av(x)  is  the  unstrained  average  valence  band,  av  and  b  are 
the  valence  band  deformation  potentials,2^  as  described  in 
Table  I. 
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Fig.  6.  The  strain  field  map  calculated  for  the  x-z  plane  for  dots  with  a  base  width  of  16.8  nm:  (a)  £x c,  (b)  eyy,  (c)  and  (d)  £b  =  £zz~ 
(£xx  +  fy)/ 2. 


The  carriers  confining  potentials  can  now  be  determined 
from 

Vs  =  Ec(x)  -  Ec(x  =  0), 

yW,  _  £hh(:c  _  o)  _ 

and 

V"1  =  E*(jc  =  0)  -  E^(x),  (22) 

where  x  =  0  describes  the  conditions  inside  the  GaAs  matrix 
(with  no  indium),  and  x^=0  reflects  the  properties  inside  the 
In^Gai_xAs  quantum  dot. 

The  material  parameters  used  in  our  calculations  are 
summarized  in  Table  I.  Several  different  models  have  been 
developed  to  calculate  the  electronic  structure  of  In^Gai-^As 
dots.14,25~27)  Here,  we  use  an  eight-band  k  p  description, 
with  strain  effects  included  via  deformation  potential 
theory.2^  The  eight-band  Luttinger-Kohn  parameters  are 
calculated  according  to  the  approach  used  by  Pollak28)  and 
Bir  et  al29)  The  bound  states  of  the  quantum  dots  are  found 
by  numerically  solving  the  Schrodinger  equation,  which,  in 
the  effective-mass  approximation,  is  written  as 


In  the  expression  above,  m*(r)  is  replaced  by  m£jGaAs(r) 
inside  the  dot,  and  by  m^As(r)  in  the  GaAs  matrix  material. 
The  parameter  V(r)  is  the  three-dimensional  confining 
potential.  In  the  frame  of  the  eight-band  model,  the  wave 
function  can  be  expanded  as 
8 

^(r)  =  ^Fv(r)«;(r).  (24) 

j=  1 

As  described  by  Stier  et  al2A)  the  k  p  model  has  a  few 
drawbacks  when  it  is  applied  to  nanostructures.  These 
drawbacks  are  related  to  the  fixed  number  of  Bloch  functions 
used  in  synthesizing  the  wavefunction  (eight  in  the  eight- 
band  formalism).  The  problems  are  also  connected  with  the 
restriction  to  the  Brillouin  zone  center,  T,  and  the 
assumption  of  the  same  Bloch  functions,  regardless  of 
material  and  strain  variations.  The  choice  of  the  appropriate 
boundary  conditions  and  the  arbitrariness  of  the  matching 
conditions  for  the  envelope  functions  across  the  heterointer¬ 
face  are  also  issues  to  be  concerned  with.  In  our  case,  we  use 
the  guidelines  developed  by  Burt30)  for  the  exact  envelope- 
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Fig.  7.  The  strain  field  map  calculated  in  the  x-z  plane  for  dots  with  a  base  width  of  19.7  nm:  (a)  (b)  %,  (c)  e and  (d)  £b  = 

(£xx  +  %)/ 2. 


Table  I.  Material  parameters  used  in  theoretical  calculations. 


Parameter 

GaAs 

InAs 

hvGai-jAs 

o(A) 

5.6503 

6.0553 

5.6503  +  0.405* 

E&  (eV) 

6.4  K 

1.518 

0.413 

1 .5 1 8  —  1 .580*  +  0.475a2 

300  K 

1.424 

0.324 

0.324  +  0.7(1  -  x)  +  0.4(1  -  *)2 

Ao  (eV) 

0.340 

0.380 

0.340  -0.093x-H0.133x2 

C„  (N/m2) 

11.88 

8.33 

11.88  -  3.55a: 

Ca  (N/m2) 

5.38 

4.53 

5.38  -  0.85* 

flc(eV) 

-8.013 

-5.08 

-8.013  +  2.933* 

av  (eV) 

1.16 

1.00 

MeV) 

-1.7 

-1.8 

(eV) 

-6.92 

-6.747 

—6.92  +  0.231*  —  0.058*2 

< 

0.0667  mQ 

0.02226  mo 

0.0667  -  0.0419*  -  0.00254*2 

function  theory  by  taking  into  account  the  variation  of  mass 
parameters  from  their  bulk  values  due  to  strain-induced  band 
deformations. 

Equation  (23)  can  be  solved  for  the  eigenvalues  and 
eigenfunctions  of  the  system  by  invoking  periodic  boundary 


Small  dots  (b=16.8  nm)  Big  dots  (b=19.7  nm) 
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Fig.  8.  Schematic  representation  of  calculated  energy  levels  (referenced 
to  the  vacuum  level)  for  quantum  dots  with  base  widths  of  16.8  nm  and 
19.7  nm,  respectively.  The  fundamental  interband  transition  energies, 
(Eco  —  £ho),  and  the  intraband  transition  energies,  (Ec\  —  Et o),  are  also 
shown. 


conditions,  expanding  ^n(r)  in  terms  of  normalized  plane- 
wave  states,  and  diagonalizing  the  resulting  matrix.  This 
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approach  has  also  been  used  by  Cusack  et  al  ;25)  it  does  not 
require  the  matching  of  the  wavefunctions  across  the 
boundary  between  the  dot  and  the  matrix  material.  This 
makes  the  method  useful  when  arbitrary  (position  depen¬ 
dent)  confining  potentials  are  involved.  The  boundary 
conditions  would  only  require  that  the  states  do  not  overlap 
significantly  for  neighboring  dots. 

The  energy  levels  of  the  truncated  pyramidaly-shaped 
Ino.23Gao.77 As  quantum  dots  have  been  calculated  using  the 
parameters  listed  in  Table  I.  The  conduction  ground  state 
energy  level  for  dots  with  a  base  of  b  =  19.7  nm  is 
calculated  to  be  Er n  =  “5.647  eV  (measured  from  the 
vacuum  level);  the  first  excited  state  for  electrons  is  located 
at  Eei  =  —5.540  eV.  For  the  heavy  holes,  the  ground  state 
energy  (measured  from  the  vacuum  level)  is  located  at 
£h0  =  —6.773  eV.  The  fundamental  transition  energy, 
(Eeo  —  Eho),  is  therefore  1.126  eV.  The  separation  between 
the  conduction  band  ground  state  and  the  first  excited  state 
has  been  calculated  as  107  meV. 

For  small  dots,  characterized  by  a  base  dimension  b  — 
16.8  nm,  the  ground  state  for  electrons  is  calculated  as  Ee0  = 
“5.474 eV  (as  measured  from  the  vacuum  level);  the  first 
excited  state  is  located  at  Et  1  =  —5.609  eV.  For  the  heavy- 
holes,  the  ground  state  energy  (measured  from  the  vacuum 
level)  is  Zq,o  =  —6.775  eV.  The  fundamental  transition 
energy,  (£co  —  EhoX  *s  therefore  1.166  eV.  The  separation 
between  the  conduction  ground  state  and  the  first  excited 
state  is  evaluated  as  135  meV.  The  conduction  ground  state 
of  the  small-size  dots  is  calculated  to  lie  about  38.4  meV 
above  the  ground  state  of  large-size  dots.  All  the  calculated 
energy  levels  are  shown  in  Table  II;  they  are  also  displayed 
in  Fig.  6  in  a  simple  energy  band  diagram.  The  results  of  the 
calculations  performed  here  are  in  good  agreement  with  the 
photoluminescence,  as  well  as  infrared  absorption  data 
obtained  at  300  K. 

We  want  to  emphasize  that  the  dot-sizes  chosen  for  the 
calculations  are  not  unique.  The  strain-modified  confined 
energy  levels  are  not  as  sensitive  to  variation  of  the  height  of 
the  dot  as  they  are  to  the  change  of  the  base  length  of  the 
structure.  Although  other  combinations  of  the  base  and 
height  might  yield  similar  ground  state  energy  levels,  a 
unique  spectrum  (ground  state  and  excited  states  for  both 
electrons  and  heavy  holes)  would  be  obtained  only  for  a 
particular  combination  of  these  dimensions.  Our  choice  of 
dimensions  in  this  work  was  guided  by  measurements 
obtained  from  transmission  electron  microscope  images  of 
the  dots.  The  final  choice  was  guided  by  the  experimental 
transition  energies  obtained  from  the  photoluminescence  and 
absorption  measurements  we  performed. 


4.  Summary 

In  summary,  we  have  presented  results  of  our  optical 
study  of  the  thermal  redistribution  of  carriers  in  quantum-dot 
ensembles  of  various  sizes.  Using  an  analytical  approach,  we 
have  calculated  the  energy  levels  in  the  dots  and  found  that 
the  theoretical  results  are  in  agreement  with  our  experi¬ 
mental  data. 

Acknowledgments 

This  work  is  supported  by  the  US  Army  Research  Office, 
Research  Triangle  Park.  North  Carolina  and  by  the  US  Army 
Research  Laboratory,  Adelphi,  Maryland. 

1)  F.  Heinrichsdorff,  Ch.  Ribbat,  M.  Grundmann  and  D.  Bimberg:  Appl. 
Phys.  Lett.  76  (2000)  556. 

2)  K.  Imamura,  Y.  Sugiyama,  Y.  Nakata,  S.  Muto  and  N.  Yokoyama:  Jpn. 
J.  Appl.  Phys.  34  (1995)  L1445  . 

3)  D.  Pan,  E.  Towe  and  S.  Kennerly:  Appl.  Phys.  Lett.  73  (1998)  1937  . 

4)  R.  L.  Sellin,  Ch.  Ribbat,  M.  Grundmann,  N.  N.  Ledentsov  and  D. 
Bimberg:  Appl.  Phys.  Lett.  78  (2001)  1207. 

5)  S.  Maimon,  E.  Finkman,  G.  Bahir,  S.  E.  Schacham,  J.  M.  Garcia  and 
P.  M.  Petroff:  Appl.  Phys.  Lett.  73  (1998)  2003. 

6)  D.  Pan,  E.  Towe  and  S.  Kennerly:  Appl.  Phys.  Lett.  76  (2000)  3301. 

7)  D.  I.  Lubyshev,  P.  P.  Gonzalez-Borrero,  E.  Marega,  Jr.,  E.  Petiprez, 
N.  la  Scala,  Jr.  and  P.  Basmaji:  Appl.  Phys.  Lett.  68  (1996)  205. 

8)  R.  P.  Mirin,  J.  P.  Ibbetson,  K.  Nishi,  A.  C.  Gosard  and  J.  E.  Bowers: 
Appl.  Phys.  Lett.  67  (1995)  3795. 

9)  S.  Fafard,  S.  Raymond,  G.  Wang,  R.  Leon,  D.  Leonard,  S. 
Charbonneau,  J.  L.  Merz,  P.  M.  Petroff  and  J.  E.  Bowers:  Surf.  Sci. 
361  (1996)  778. 

10)  L.  Chu,  A.  Zrenner,  G.  Bohm  and  G.  Abstreiter:  Appl.  Phys.  Lett.  75 
(1999)  3599. 

11)  J.  Phillips,  P.  Bhattacharya,  S.  W.  Kennerly,  D.  W.  Beckman  and  M. 
Datta:  IEEE  J.  Quantum  Electron.  35  (1999)  936. 

12)  A.  Weber,  O.  Gauthier-Lafaye,  F.  H.  Julien,  J.  Brault,  M.  Gendry,  Y. 
Desieres  and  T.  Benyatton:  Appl.  Phys.  Lett.  74  (1999)  413. 

13)  P.  Hawrylak,  A.  Wojs,  D.  J.  Lockwood  and  P.  D.  Wang:  Surf.  Sci.  361 
(1996)  774. 

14)  M.  Grundmann,  O.  Stier  and  D.  Bimberg:  Phys.  Rev.  B  52  (1995) 
11969. 

15)  L.  R.  C.  Fonseca,  J.  L.  Jimenez  and  J.  P.  Leburton:  Phys.  Rev.  B  58 
(1998)  9955. 

16)  H.  Xu,  Q.  Gong,  B.  Xu,  W.  Jiang,  J.  Wang,  W.  Zhou  and  Z.  Wang:  J. 
Cryst.  Growth  200  (1999)  70. 

17)  Y.  P.  Varshini:  Physica  34  (1967)  149. 

18)  A.  Polimeni,  A.  Patane,  M.  Henini,  L.  Eaves  and  P.  C.  Main:  Phys. 
Rev.  B  59  (1999)  5064. 

19)  A.  Patane,  A.  Polimeni,  P.  C.  Main,  M.  Henini  and  L.  Eaves:  Appl. 
Phys.  Lett.  75  (1999)  814. 

20)  V.  G.  Stoleru,  D.  Pal  and  E.  Towe:  submitted  to  Physica  E  (2001). 

21)  J.  Eshelby:  Proc.  R.  Soc.  London,  Ser.  A  241  (1957)  376. 

22)  J.  R.  Downes,  D.  A.  Faux  and  E.  P.  O’Reilly:  J.  Appl.  Phys.  81  (1997) 
6700. 

23)  G.  Cipriani,  M.  Rosa-Clot  and  S.  Taddei:  Phys.  Rev.  B  61  (2000) 
7536. 


Table  II.  Theoretically  calculated  energy  levels,  along  with  the  experi¬ 
mentally  determined  interband  emission  and  intraband  absorption 
energies. 

Base  width  =  19.7  nm  Base  width  =  16.8  nm 
Interband  Intraband  Interband  Intraband 

emission  absorption  emission  absorption 

Experimental  1.12eV  107meV  1.16eV  135meV 

Theoretical  1.126eV  107meV  1.166eV  135meV 


24)  O.  Stier,  M.  Grundmann  and  D.  Bimberg:  Phys.  Rev.  B  59  (1999) 
5688. 

25)  M.  Cusack,  P.  R.  Briddon  and  M.  Jaros:  Phys.  Rev.  B  54  (1996) 
R2300. 

26)  J.  Kim,  L.-W.  Wang  and  A.  Zunger:  Phys.  Rev.  B  57  (1998)  R9408. 

27)  J.-Y.  Marzin  and  G.  Bastard:  Solid  State  Commun.  92  (1994)  437. 

28)  F.  H.  Poliak:  Semicond.  Semimet.  32  (1990)  17. 

29)  G.  L.  Bir  and  G.  E.  Pikus:  Symmetry  and  Strain-Induced  Effects  in 
Semiconductors  (Wiley,  New  York,  1974). 

30)  M.  G.  Burt:  J.  Phys.:  Condens.  Matter  4  (1992)  6651. 


